AD-H147  311  A  REVIEW  OF  IMPERFECTION  SENSITIVITY  OF  STIFFENED  1/2 

SHELLSCU)  ANAMET  LABS  INC  SAN  CARLOS  CA  APPLIED 
MECHANICS  DIV  R  L  CITERLEY  FEB  84  ANAMET-ASIAC-1183.  1R 
UNCLASSIFIED  AFHAL-TR-84-3B08  F33615-81-C-3201  F/G  13/13  NL 


I1, 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


Af  WAL-TR-84-3006 


A  REVIEW  OF  IMPERFECTION  SENSITIVITY  OF  STIFFENED  SHELLS 


Anamet  Laboratories,  Inc. 
Applied  Mechanics  Division 
San  Carlos,  California 


February  1984 


■r  i  ll  Report  for  Period  December  1982  -  November  1983 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


Fliqht  Dynamics  Laboratory 

Air  Force  Wright  Aeronautical  Laboratories 

Air  Force  Systems  Command 

Wright-Patterson  Air  Force  Base,  Ohio  45433 


s 


ELECTE 
NOV  9  1984 


3 


/  B 


10  2  3  1*?® 


reproduced  at  government  expense 


NOTICE 

When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  related  Government  procure¬ 
ment  operation,  the  United  States  Government  thereby  incurs  no  responsibility 
nor  any  obligation  whatsoever;  and  the  fact  that  the  government  may  have  formu¬ 
lated,  furnished,  or  in  any  way  supplied  the  said  drawings,  specifications,  or 
other  data,  is  not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner 
licensing  the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in 
any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


Project  Engineer 


Chief,  Analysis  &  Optimization  Branch 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  your  organization,  please  notify 

AFMAL/FIBRA _ _ ,  Wright-Patterson  AFB,  OH  45433,  to  help 

us  maintain  a  current  mail ing  list. 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obi Jons,  or  notice  on  a  specific  document. 


UNCLASSIFIED _ 

s£CUBlT\  classification  OF  this  page 


la  REPORT  SECURITY  CLASSIFICATION 


luitfvivtaiaiaii 


2a  SECURITY  CLASSIF  i  CATION  AUTHORITY 


20  DECLASSIFICATION  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBERlS) 

AS  I AC  REPORT  NO.  1183. 1A 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL 

(If  applicable) 

ANAMET  LABORATORIES,  INC. 


6c  AOORESS  .Cify.  State  and  ZAP  Code ) 

100  INDUSTRIAL  WAY 

SAN  CARLOS,  CALIFORNIA  94070 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


S«  NAME  OF  F  UNDlNG-'SPONSORlNG 
ORGANIZATION 


8c  ADDRESS  -City.  State  and  ZIP  Code ) 


8b.  OFFICE  SYMBOL 
(If  applicable  i 


11  TITLE  -lnctude  Security  Classification) 

A  REVIEW  OF  IMPERFECTION  SENSITIVITY  OF 


12  personal  authorisi  STIFFENED  SHELLS  (U 
Citerlev,  Richard  L. 


13b.  time  COVERED 


3.  distribution/availability  of  report 

APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


5.  monitoring  organization  report  numberisi 

AFWAL-TR-84-3006 


7a.  NAME  OF  MONITORING  ORGANIZATION 

FLIGHT  DYNAMICS  LABORATORY  (AFWAL/FIBRA) 


7b.  AOORESS  (City.  State  and  ZIP  Code) 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
WRIGHT-PATTERSON  AFB,  OHIO  45433 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F33615-81-C-3201 


10.  SOURCE  OF  FUNDING  NOS 


PROGRAM 
E  LEMENT  NO. 

62201F 


PROJECT 

TASK 

WORK  UNIT 

NO. 

NO. 

NO. 

2401 

02 

45 

13a.  T>  PE  OF  REPORT 


14.  DATE  OF  REPORT  (Yr.  Mo..  Day I  I  15  PAGE  COUNT 


ii'.wi*:! 


1984  February 


COS AT  I  CODES 


18  SUBJECT  TERMS  < Continue  on  reverte  if  neceuary  and  identify  by  block  numbert 

Shells  (Structural  Forms) 

Bifurcation 

Static  Stability 

Imperfection  Sensitivity 

Buckl ina 

_ Stiffened  Shells _ 

19  A8STRACT  Continue  on  revert e  if  neceuar y  and  identify  by  block  number i 

‘A  short  history  is  related  about  how  imperfections  (i.e.,  deviations  from  the  assumed 
theoretical  form  of  a  shell  structure)  came  to  be  regarded  as  the  major  contributor  to 
the  experimental  scatter  observed  in  shell  buckling  tests,  and  the  discrepancy  between 
experimentally  determined  buckling  loads  and  those  determined  by  classical  linear  theory. 
J.M.T.  Thompson's  insight  into  the  relationship  between  imperfection  sensitivity  and 
optimality  of  shell  structures  with  respect  to  load  configurations  is  brought  out. 

The  significant  contributions  to  understanding  elastic  postbuckling  behavior  of  thin 
stiffened  shells  and  how  this  work  relates  to  imperfection  sensitivity  are  reviewed. 

The  implications  of  this  work  and  its  applicability  to  current  engineering  design  and 
analysis  are  assessed.  The  effects  of  geometric  imperfections  as  well  as  those  induced 
by  boundary  conditions  or  multiple  load  cases  on  the  postbuckling  behavior  of  shells  are 
discussed.  Added  discussions  by  Profs.  Singer  and  van  Neut  are  also  presented. 

I 


20  OISTRIBUTION/AVAILABILITY  of  abstract 
UNCLASSIFIED  UNLIMITED  QC  SAME  AS  RPT  □  OTIC  USERS  □ 


22*  NAME  OF  RESPONSIBLE  INOIVIOUAL 

■Iannis  R.  Johnson 


00  FORM  1473,  83  APR  EDITION  OF  1  JAN 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b  TELEPHONE  NUMBER 
(Include  Area  Code) 

513-2556992 


22c  OFFICE  SYMBOL 

AFWAL/F1BR 


EDITION  OF  1  JAN  73  IS  OBSOLETE. 


UNCLASSIED _ 

SECURITY  CLASSIFICATION  ■••F  THIS  PAGE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  this  page 


BLOCK  18  (continued) 

Postbuckl ing 
Shell  Buckling 
Buckling  Loads 


Elastic  Stability 
Structural  Stability 
Unstiffened  Shells 
Catastrophe  Theory 


TABLE  OF  CONTENTS 


Page 


1.0  INTRODUCTION  .  1 

2.0  FUNDAMENTAL  PRINCIPLES  AND  THEIR  HISTORICAL 

DEVELOPMENT  .  4 

2.1  Basic  Concepts  . .  4 

2.2  Semi -Empirical  Methods  . . 17 

2.3  Geometric  Nonlinear  Theories  .  21 

3.0  RECENT  DEVELOPMENTS  IN  ANALYTICAL  METHODS  .  22 

3.1  Unstiffened  Shells  .  23 

3.1.1  Equilibrium  Method  .  23 

3.1.2  Koiter's  Method  .  27 

3.1.3  Thompson's  Extension  of 

Koiter's  Method  .  32 

3.2  Imperfection  Sensitivity  of 

Stiffened  Shells  .  40 

3.2.1  Equilibrium  Methods  .  41 

3.2.2  Koiter's  Energy  Methods  .  48 

4.0  NUMERICAL  METHODS  .  58 

5.0  EXPERIMENTS  .  60 

6.0  CONCLUSIONS  AND  SUMMARY  .  67 

REFERENCES  .  71 

DISCUSSION  .  99 

AUTHOR'S  CLOSURE .  102 

ADDITIONAL  REFERENCES  .  105 


iii 


y*. 


u 


LIST  OF  ILLUSTRATIONS 


Figure  No. 


Page 


Load  deflection  curve-effects 
of  imperfections  . . 


Typical  postbuckling  paths  for  cylinders  ....  12 

Response  of  cylinder  to  axial  load  .  13 

Equilibrium  paths  with  stability  loss 

at  branch  points  .  1 4 

Possible  bifurcation  paths  . .  16 

Two-mode  branching  . 25 

Effect  of  axisymmetric  imperfection  on 
buckling  of  cylindrical  shell  .  30 


Uncoupled  branching  configurations  for 
doubly-symmetric  structural  systems  .  36 

Forms  of  coupled  postbuckling  for  ideal 
doubly-symmetric  structural  systems  .  37 

Forms  of  postbuckling  equilibrium  path 

for  doubly-symmetric  structural  systems 

with  imperfections  .  37 

Bifurcation  paths  based  upon 

catastrophe  theory  .  39 

Typical  postbuckling  behavior  of  axially 

loaded  orthotropic  shell  .  43 

Postbuckling  functions  of  stiffened 

cylinders  axially  loaded  .  44 


Postbuckling  functions  of  axially 

compressed  pressurized  orthotropic 

cylinders  -  axially  stiffened  .  45 

Postbuckling  functions  of  axially 

compressed  pressurized  orthotropic 

cylinders  -  circumferentially  stiffened  .  46 


v 


LIST  OP  ILLUSTRATIONS  (Concluded) 

Fictile  No.  Page 

16  Classical  buckling  and  imperfection- 
sensitivity  of  simply  supported,  axially 
stiffened  cylinders  under  axial 

compression  . 52 

17  The  effect  of  stringer  eccentricity  on 

the  buckling  and  postbuckling  behavior 
of  axially  stiffened  cylindrical  shells 
which  are  simply  supported  at  the  skin 
middle  surface  and  loaded  in  axial 
compression  . 55 

18  The  effect  of  stringer  eccentricity  on  the 
buckling  and  postbuckling  behavior  of  clamped 
axially  stiffened  cylindrical  shells  which 


are  loaded  in  axial  compression  .  56 

19  Early  experimental  results  for 

clamped  shells  .  6l 

20  Summary  of  sphere  test  performance  .  63 

A-l  Comparative  surfaces  for  catastrophe 

theory  and  energy  methods  of  a  shallow 
spherical  cap  .  10A 


LIST  OP  TABLES 


1.0  INTRODUCTION 


The  stiffened  shell  is  perhaps  the  most  utilized  structure 
in  industry  today.  Since  the  early  1930’s,  the  aerospace  and 
aircraft  industry  has  used  this  class  of  shell  to  its  fullest 
extent.  However,  as  the  shell  structure  is  made  thinner  and 
lighter  for  a  given  loading,  it  becomes  prone  to  a  catastrophic 
collapse  known  as  shell  instability.  The  magnitude  of  the  load 
that  causes  the  shell  to  become  unstable  is  significantly 
affected  by  the  size  and  geometry  of  imperfections  in  the  initial 
shape  of  the  shell.  Consequently,  in  order  to  ensure  an  effi¬ 
cient  design,  considerable  research  has  been  conducted  in  the 
United  States  and  overseas  to  more  fully  understand  the  complex 
behavior  of  these  structures. 

The  development  of  an  understanding  of  the  behavior  of  shell 

structures  had  its  start  around  the  turn  of  the  century.  Surveys 

by  Flugge10^  ,  Langhaar1^,  Nowinski22\  Naghdl21^  and 
ppn  ppi 

Nash  *  clearly  demonstrate  that  even  after  more  than  three 
quarters  of  a  century  of  study,  there  is  a  continued  interest  in 
shell  analysis  and,  in  particular,  elastic  shell  stability 
analysis.  In  fact,  the  interest  has  been  so  great  that  the  total 
number  of  papers  written  on  the  subject  exceeds  six  thousand,  of 
which  nearly  one  thousand  were  written  in  the  last  decade. 

The  interest  in  shell  analysis  has  not  only  been  associated 
with  the  aerospace  industry,  but  also  in  the  marine  and  energy 
fields,  just  to  mention  a  few.  For  example,  the  Subcommittee  on 
Shells  for  the  Pressure  Vessel  Research  Committee  of  the  Welding 
Research  Council  has  supported  several  studies  for  the  analysis 
of  shell  structures.  Also,  a  chronicle  of  support  of  the  Pres¬ 
sure  Vessel  and  Piping  Division  of  the  ASME  has  been  collected  in 
several  ASME  documents'^* The  author  has  also  presented  a 
limited  survey  of  the  imperfection  sensitivity  and  postbuckling 


•Superscript  numbers  denote  references  found  at  end  of 
report . 
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of  shells  as  part  of  the  Pressure  Vessel  and  Piping  Decade  of 
Progress  of  1982?9.  a  more  in-depth  review  of  buckling  of  shells 
as  a  general  category  has  been  given  by  Bushnell^ » . 

After  all  that  has  been  written  on  the  subject,  one  might 
assume  that  the  subject  matter  has  been  amply  covered.  Some 
researchers  have  taken  this  attitude,  but  many  design/analyst 
engineers  find  that  they  must  still  rely  upon  ingenuity  and  engi¬ 
neering  instincts  when  evaluating  the  stability  characteristics 
of  shell  structures  which  are  subject  to  a  specific  complex  load¬ 
ing  condition.  Perhaps  the  present  state-of-the-art  is  too  math¬ 
ematical  or  abstract  for  the  practicing  engineer.  Perhaps  the 
methods  of  analysis  are  incomplete  with  respect  to  the  under¬ 
standing  of  the  instability  of  shells.  This  is  more  eloquently 
presented  by  Kolter1^  in  his  discussion  of  Truesdell  and  Noll's 
remarks^0*}  on  the  purpose  of  the  elastic  stability  literature. 
Koiter  states: 

"Our  first  aim,  to  show  the  real  purpose  of  existing  theo¬ 
ries  and  current  research  in  elastic  stability,  is  by  far  our 
easiest  task.  Structural  stability  is  in  fact  one  of  the  most 
important  criteria  in  the  design  of  many  engineering  structures. 
Our  knowledge  in  this  field  is  admittedly  far  from  complete,  in 
particular  for  structures  in  which  inelastic  behaviour  or  dynamic 
loading  conditions  are  essential  features,  but  the  theory  of 
elastic  stability,  however  incomplete  it  may  be,  constitutes  an 
indispensable  tool  in  achieving  properly  designed  and  efficient 
structures,  capable  of  withstanding  their  service  loads  without 
catastrophic  failures.  It  is  quite  true,  of  course,  that  many 
investigations  of  elastic  stability  lean  heavily  on  more  or  less 
crude  approximations  and  on  heuristic  physical  arguments,  not  by 
choice,  however,  but  only  due  to  the  absence  of  a  rigorous  theory 
which  is  at  the  same  time  capable  of  a  reasonably  simple  applica¬ 
tion  to  the  problem  at  hand.  Elastic  stability  shares  the  fate 
of  all  engineering  science  that  it  cannot  afford  to  await  the 
development  of  fully  rigorous  mathematical  theories  before  it 
deigns  to  consider  the  solution  of  stability  problems  arising  in 


_ _ 
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engineering  practice.  Its  servant’s  role  in  engineering  implies 
that  it  must  try  and  solve  these  problems,  if  possible  by  rigor¬ 
ous  analysis,  if  necessary  by  hook  or  by  crook.  A  rough  estimate 
today  may  be  worth  more  than  an  accurate  approximation  a  month 
from  now." 

At  present,  the  dilemma  that  the  practicing  engineer  faces 
is  that  he  has  only  two  choices  in  examining  the  stability  of 
imperfect  shells  —  he  must  use  either  (1)  complex  computer  pro¬ 
grams  which  are  based  upon  certain  theories  and  therefore  have 
certain  limitations  —  some  specified,  others  not;  or  (2)  simpli¬ 
fied  equations  modified  by  semi-empirical  (graphical)  data 
developed  from  experiments  that  may  or  may  not  have  been  ade¬ 
quately  controlled.  Without  a  complete  understanding  of  the 
subject  matter,  the  engineer  could  easily  misuse  these  computer 
programs  and/or  data,  and  misleading  conclusions  could  be  drawn. 
The  dilemma  is  further  aggravated  because  modern  technology 
demands  that  more  efficient  structures  be  designed,  sometimes 
beyond  the  imposed  limitations  of  presently  accepted  data  or 
computational  techniques. 

It  is  therefore  the  purpose  of  this  review  to  illustrate  the 
basis  of  the  accepted  theories  of  the  imperfection  sensitivity  of 
shells  and  also  to  suggest  works  that  may  be  helpful  to  the 
practicing  engineer  in  gaining  additional  insight.  Since  the 
previously  mentioned  surveys  give  a  rather  broad  description  and 
chronicle  of  shell  instability  studies.  It  would  appear  to  be 
appropriate  for  this  review  to  present  material  which  permits  the 
practicing  engineer  to  focus  on  the  complexities  involved  in  the 
so-called  elementary  stability  theories  and  understand  their 
Limitations.  In  this  way  he  may  bridge  the  gap  between  his  needs 
and  those  accepted  theories.  Because  of  the  enormous  volume  of 
papers  on  just  one  subject,  it  may  require  the  engineer  to  review 
several  papers  to  gain  the  required  understanding  of  a  single 
concept.  A  summary  Is  provided  to  aid  in  this  evaluation.  The 
review  must  necessarily  present  abbreviated  forms  of  some  rather 
cumbersome  mathematical  expressions.  Their  appearance  is  by  no 


means  intended  to  be  complete,  but  rather  a  reminder  of  their 
complexity. 

2.0  FUNDAMENTAL  PRINCIPLES  AND  THEIR  HISTORICAL  DEVELOPMENT 
2.1  Basic  Concepts 

The  foundation  of  the  analysis  of  thin  shells  of  revolution 
was  developed  more  than  a  century  ago  and  Novozhilov22^  presented 
a  rather  comprehensive  description  of  this  development.  The 
basic  equations  of  equilibrium  were  developed  by  Aron2^  in  187^ 
and  Love-*-95  in  1888.  In  these  early  developments  it  was  shown 
that  the  main  feature  that  distinguishes  thin  shell  theory  from 
the  three-dimensional  theory  of  elasticity  is  that  the  shell  can 
be  regarded  as  a  two-dimensional  surface.  All  fundamental  vari¬ 
ables  are  therefore  dependent  on  two  orthogonal  curvilinear 
coordinates,  and  c2,  which  define  a  point  on  a  middle  or 
reference  surface.  The  description  of  this  reference  surface 
obeys  the  relationships  of  differential  geometry.  The  governing 
equations  can  be  greatly  simplified  if  the  reference  surface  is 
restricted  to  a  surface  of  revolution.  The  key  parameters  that 
describe  the  shell  of  revolution’s  geometry  are  its  two  principal 
radii  of  curvature,  and  R2.  For  practical  shell  applications, 
the  coordinate  system  usua-'-^y  oriented  along  the  princi¬ 

pal  curvature  directions.  For  a  shell  of  revolution,  the  in¬ 
plane  coordinates  are  in  the  meridional,  where  d^  =  R^d<|>,  and 
the  circumferential,  d£^  =  Rgde,  ^^rec^^ons*  The  third  ortho¬ 
gonal  coordinate  is  direction  and  is  in  the  direction  of  the 
outward  normal. 

The  development  of  the  strain-displacement  relations  for  a 
general  shell  in  terms  of  the  orthogonal  curvilinear  coordinates 
is  given  by  Love.  These  equations,  as  well  as  the  remaining 

1  ft  P 

fundamental  equations  of  thin  shells,  are  presented  by  Krausx  . 
in  his  presentation,  the  accepted  definitions  are  used  to  illus¬ 
trate  Reissner's  version  of  the  Love  theory.  Four  representive 
shell  theories  are  discussed  which  preserve  Love’s  original 
assumptions . 


The  measure  of  deformation  of  the  reference  surface  can  be 
expressed  through  a  set  of  strains  and  curvatures  which  are  given 
in  terms  of  the  displacement  vector  of  the  middle  surface  and  the 
geometric  parameters  of  the  surface.  The  reference  surface 
strain  and  curvature  create  stresses  throughout  the  shell.  Inte¬ 
grating  the  stresses  and  their  first  moments  about  the  reference 
surface  give  the  in-plane  or  membrane  stress  resultants  and  the 
bending  stress  resultants.  Consistent  with  the  theory  of  Love's 
first  approximation  are  the  following  four  assumptions: 

(1)  The  thickness  of  the  shell,  t,  is  small  such  that 
t/R^  and  t/R2  <<  1. 

(2)  The  strains  and  displacements  are  small  so  that  second 
order  terms  and  higher  may  be  neglected. 

(3)  The  extensional  stresses  normal  to  the  middle  surface 
may  be  considered  small  with  respect  to  the  other 
stress  components. 

(4)  Normals  to  the  reference  surface  prior  to  deformation 
remain  normal  after  deformation  and  have  no  extension. 

The  displacement  of  a  point  on  the  reference  surface  of  the 
shell  Is  expressed  by  u,  v  and  w  in  the  and  directions, 

respectively;  and  the  displacement  of  a  point  of  a  distance,  Z, 
away  from  the  reference  surface  can  be  given  by: 

A 

ui  =  u  +  ZfJi 


u  =  v  +  Z  B 

U  ' 


w  =  w 


where  b  and  B_  are  the  reference  rotation  about  the  surface 

<j>  y 

coordinates  anfi  respectively.  Some  theories  permit  local 

rotation  (twist,*)  about  the  normal  coordinate  as  well. 

In  order  to  fully  understand  the  implication  of  Love's  first 
approximation,  the  equations  of  equilibrium  and  the  strain- 
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displacement  equations  must  be  considered  in  conjunction  with  the 
relationship  between  the  stresses  and  strains,  such  as  Hooke's 
Law  for  a  linear  elastic  material.  Love's  first  assumption 
defines  what  is  meant  by  thin.  Although  a  precise  definition  is 
not  available,  usually  a  radius  to  thickness  ratio  R/t  >  10  has 
been  accepted  by  many.  The  second  assumption  assures  linearity 
of  the  governing  differential  equations.  The  third  and  fourth 
assumptions  require  the  normal  stress  and  strain  and  the 
transverse  shear  stresses  to  be  negligible.  Usually,  for  homo¬ 
geneous,  orthotropic  materials,  there  are  three  mutually  per¬ 
pendicular  planes  of  elastic  symmetry,  and  these  planes  also 
usually  lie  in  the  planes  of  principal  curvature.  Thus,  the 
theories  are  applicable  to  orthotropic  shells  of  revolution. 

Hildebrand,  Reissner  and  Thomas^^  expand  on  this  theme  and 
give  a  more  complete  description  of  the  higher  order  approximate 
theories.  Several  other  methods  have  been  presented  by  a  number 
of  investigators  in  an  attempt  to  improve  the  analysis  of  shells. 
An  excellent  discussion  of  the  various  derivations  of  the  govern¬ 
ing  equations,  their  similarities  and  their  differences,  is  given 
by  Koga  and  Endo1^,  and  by  Leissa1^2.  Over  fourteen  different 
theories  are  examined  by  Leissa.  These  include  Novozhllov,  Love, 
Timoshenko^O^,  Byrne^,  Flugge^^,  Goldenveizer -*-20,  Lur'ye^^ 
Keissner2^0,  Vlasov^^  Sanders2^,  Donnell^,  and  Mushtarl2^ . 

In  Leissa 's  comparison,  he  shows  that  the  usual  result  when 
adopting  any  of  these  theories  is  that  inconsistent 
approximations  are  required  in  order  to  make  any  significant 
improvements  over  Love's  first  approximation. 

With  the  fundamental  linear  shell  equations  addressed,  at¬ 
tempt.?  were  made  between  1900  and  1914  by  Lorenz ^94, 
von  Mlses^1',  Southwell2®^-  and  Timoshenko^*^  to  examine  the  crit¬ 
ical  buckling  load  of  a  thin  circular  cylindrical  shell  under 
axial  compression.  The  linear  equations  governing  the  buckling 
of  a  cylinder  in  a  membrane  prestressed  state  of  compression  due 
to  an  axial  load,  P,  can  be  given  in  the  form: 


E1(u,v,w)  =  0 


(2) 


E2(u,v,w)  =  0 
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where  E^  and  E2  are  equilibrium  equations  in  the  and  c2  direc¬ 
tions  in  terms  of  u,  v,  and  w;  R  is  the  radius  of  the  cylinder,  x 
is  the  meridional  coordinate  (synonymous  with  the  <p  direction), 
a  is  the  circumferential  independent  variable,  and  commas  denote 
partial  differentiation.  The  variables  C  and  D  are  given  by: 
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where  E  is  the  modulus  of  elasticity,  v  is  Poisson's  ratio.  Also 
the  biharmonic  operator  is  defined  as: 
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These  three  governing  equations  described  in  Eq.  (2)  constitute  a 
linear  eigenvalue  problem,  and  the  eigenvalue  is  related  to  P. 

As  P  is  Increased  from  zero,  the  cylinder  will  undergo  an  end 
shortening  linearly  related  to  P  and  accompanied  by  a  uniform 
radial  expansion  due  to  the  Poisson  effects.  When  the  minimum 
eigenvalue  is  reached,  the  cylinder  displacement  pattern  will  be 
at  a  bifurcation  point.  For  larger  values  of  P,  the  shell  could 
continue  the  same  deformation  pattern  as  before,  or  it  could 
exhibit  a  non-uniform  radial  buckling  displacement  that  could  be 
quite  large.  The  linear  eigenvalue  formulation  does  not  allow 
any  study  of  the  post-buckling  behavior.  The  only  information  it 
provides  is  for  minimum  bifurcation  buckling  load  and  the 
associated  buckling  mode  shape. 

When  experimental  results  became  available  it  was  noted 
that,  although  good  agreement  could  be  attained  for  structures 


comprised  of  columns  and  plates,  poor  correlations  between  the 
analytical  and  experimental  buckling  loads  were  obtained  for  the 
thin-walled  circular  cylindrical  shell^^,  jn  attempts  to  recon¬ 
cile  the  discrepancy  between  experiment  and  theory  for  cylln- 
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drical  shells,  Donnell71  and  Plugge  suggested  certain 

boundary  condition  restraints  at  the  onset  of  buckling,  but  they 
could  not  account  for  the  differences  in  the  predicted  and 
observed  buckled  shapes.  Donnell's  treatment  in  1934^  suggested 
the  employment  of  a  large  deflection  theory,  which  included  an 
initial  geometric  imperfection  in  the  shape  of  the  cylinder.  The 
imperfection,  denoted  by  w,  changes  Eq.  (2)  into: 
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This  new  formulation  is  not  an  eigenvalue  problem,  and  hence  the 
cylinder  will  begin  to  deform  in  a  general  way,  with  bending,  as 
soon  as  P  is  increased  from  zero.  The  imperfection  creates  a 
pseudo  load.  Von  Karman  and  Tsien^11*  >  316  demonstrated  that 
curved  shells  can  develop  buckled  deflections  that  are  different 
from  the  original  buckling  mode  with  an  Imperfection  amplitude 
that  is  several  times  the  wall  thickness.  Von  Karman  used  the 
mode  shape  data  taken  from  the  Donnell  experiment  to  prove  his 
point  and  suggested  that  states  of  equilibrium  could  exist  at 
loads  lower  than  classical. 

During  the  early  periods  of  research  for  cylindrical  shells 
under  axial  load,  a  single  path  to  shell  collapse  (limit  or  bi¬ 
furcation  point)  was  accepted.  Postbuckllng  studies  indicated  a 
dramatic  decrease  in  the  load-carrying  ability  occurs  at  the 
bifurcation  point  and  the  cylinder  undergoes  a  dramatic  collapse. 
Figure  1  illustrates  Donnell  and  Wan's  results^®  for  an  imperfect 
cylinder  under  axial  compression.  For  many  cases  where  cylinders 
of  differing  thinness  ratios  (R/t)  were  tested,  the  normalized 
load-deflection  curve  of  the  shell  (in  terms  of  the  classical 
load  k*  and  strain  e^)  essentially  follows  a  linear  curve  until 
some  limit  point  (p)  is  reached,  and  then  the  shell  snaps  into  a 


curve-effects  of  imperfect! 


postbuckled  state.  As  the  magnitude  of  the  Imperfection  is 
Increased  (in  terms  of  an  uneveness  factor,  U),  a  marked 
reduction  in  the  maximum  load-carrying  capacity  of  the  shell  is 
experienced.  The  imperfection  can  become  so  large  that  no  snap- 
like  behavior  is  experienced  and  the  shell  simply  has  an  increas¬ 
ing  nonlinear  behavior  with  load,  as  indicated  by  the  curve  for 
UR/t  =  0.4.  The  variation  of  the  load-carrying  capacity  is  much 
less  sensitive  to  imperfections  in  the  postbuckled  regions. 
However,  the  theory  of  Donnell  and  Wan  does  predict  that  as  the 
thinness  ratio  of  the  shell  is  increased,  greater  imperfection 
sensitivity  is  experienced. 

When  a  cylindrical  shell  is  in  a  pure  membrane  state,  the 
amount  of  energy  stored  in  the  shell  can  be  quite  large  without  a 
great  deal  of  deformation.  If  the  shell  is  in  pure  axial  com¬ 
pression  prior  to  buckling,  the  shell  can  fail  (buckle)  by  devel¬ 
oping  a  rapid  exchange  of  energy  from  the  membrane  to  the  bending 
state.  This  is  usually  accompanied  by  a  large  amount  of  defor¬ 
mation.  Prebuckling  deformations  due  to  imperfections  cannot 
generally  be  seen  by  the  naked  eye,  and  yet  they  have  a  pro¬ 
nounced  effect  on  the  load-carrying  capacity  of  the  shell.  These 
deformations,  correctly  accounted  for  by  Donnell  and  Wan,  result 
from  nonlinear  considerations.  In  the  nonlinear  case,  Eq.  (4) 
takes  the  form: 
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With  the  success  of  the  predicted  behavior  of  a  cylindrical 
shell  subjected  to  an  axial  load,  other  attempts  were  made  to 
examine  the  cylindrical  shell  subjected  to  other  membrane  states 
such  as  torsion  (T)  and  hydrostatic  pressure  (p).  The  cylinder 


exhibited  imperfection  sensitivity  for  these  loading  conditions, 
but  to  a  lesser  degree  than  under  axial  load  (P),  especially  in 
the  postbuckled  region  as  shown  in  Figure  2  in  terms  of  axial 
displacement  (D),  twist,  <J»,  and  radial  displacement  WR, 
respectively. 

The  radial  displacement  satisfying  one  solution  of  the 
nonlinear  equation  can  be  expressed  in  terms  of  a  double  Fourier 
series : 

w  “  l  1  A.  .  cos  i  ~  cos  j  ^  (6) 

ij  J  x  *y 

The  number  of  arbitrary  parameters  A^j  used  in  representing  the 
postbuckled  displacement  pattern  will  dictate  what  minimum 
bifurcation  load  one  can  predict.  The  quantities  and  i  are 
the  half  wave  lengths  of  the  buckles  in  the  axial  and  circumfer¬ 
ential  directions.  The  initial  postbuckling  state  is  represented 
by  a  checkerboard  pattern  (Figure  3),  with  A  being  the  load  para¬ 
meter  and  6  being  the  axial  deflection.  When  internal  pressure 
is  considered,  the  solution  can  be  greatly  altered  such  that  the 
final  postbuckled  state  could  be  either  diamond  shaped  (for  low 
Internal  pressure)  or  ring  shaped  (for  high  internal  pressure). 
Thus,  three  possible  postbuckling  patterns  are  possible  at  the 
critical  limit  point. 

From  these  extremes  in  static  behavior  of  shells  when  sub¬ 
jected  to  a  compressive  state,  the  concept  of  a  classical  bifur¬ 
cation  buckling  can  be  developed.  At  any  load  parameter,  A, 
below  the  lowest  bifurcation  load,  the  total  energy  of  the  system 
is  increasing  for  any  small  perturbed  load.  If  generalized  co¬ 
ordinates  q^  and  qj  are  selected  such  that  they  correspond  to 
non-critical  principal  coordinates,  then  a  test  of  the  distinct 
critical  paths  can  be  made.  When  the  energy  is  increasing  with 
respect  to  any  generalized  coordinate,  the  shell  is  stable.  This 
can  be  represented  as  though  a  ball  is  resting  in  a  cup  (Point  A, 
Figure  4).  As  the  state  of  stress  approaches  infinitesimally 
closer  to  the  bifurcation  load,  a  state  of  neutral  equilibrium  is 
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reached.  Here,  the  energy  Is  conditionally  stable  to  some 
perturbed  load  and  may  be  neutrally  stable  to  a  particular 
load.  To  determine  the  stability  characteristics  of  the  shell 
beyond  this  point,  the  energy  must  be  traced  for  that  perturbed 
load.  This  action  is  as  though  the  ball  were  resting  in  a  trough 
(Point  B).  Eventually  any  further  increase  in  load  will  reach  a 
peak  (Point  C)  and  ultimate  collapse  is  reached. 

When  the  generalized  coordinate  corresponds  to  the  critical 
principal  coordinate,  postbuckling  paths  are  possible  (see 
Figures  5a-5c).  The  bifurcation  process  was  considered  to  be  a 
change  from  one  equilibrium  state  of  stress  (path  I)  to  an 
adjacent  state  of  equilibrium  (path  II).  When  the  critical 
classical  buckling  load  (denoted  as  X*)  is  reached,  infinitesimal 
adjacent  states  of  equilibrium  at  the  same  load  could  exist. 

The  nonlinear  postbuckling  response  of  a  totally  enclosed 
liner  or  oil  casing  is  represented  in  Figure  5a.  The  system  is 
stable  for  one  direction  of  motion  and  unstable  in  the  opposite, 
thus  unsymmetric  postbuckling  properties  are  possible.  The 
response  of  a  cylindrical  shell  subjected  to  axial  compression  or 
a  spherical  shell  subjected  to  external  pressure  Is  represented 
by  Figure  5b.  The  postbuckling  behavior  of  a  flat  plate  is 
represented  by  Figure  5c.  For  each  case  considered,  the  dotted 
line  represents  the  behavior  of  an  imperfect  shell. 

To  determine  the  classical  load  X*,  usually  a  set  of  ordin¬ 
ary  differential  equations  are  written  and  an  eigenvalue  of  the 
system  of  equation  is  determined.  The  governing  equations  are 
derived  from  either  equilibrium  methods,  energy  methods,  or  im¬ 
perfection  methods. 

Equilibrium  methods  can  be  used  to  examine  the  governing 
equations  for  a  structure.  Then,  adjacent  equilibrium  states  are 
hypothesized  that  satisfy  the  equations  of  equilibrium  and  the 
boundary  conditions.  These  may  be  easily  prescribed  for  simple 
systems  such  as  columns  but  may  be  impossible  to  describe  for 
complex  shell  structures. 


Figure  5.  Possible  bifurcation  paths. 


For  energy  methods,  a  deformation  state  is  assumed  to  exist 
that  satisfies  the  boundary  conditions  and  continuity.  By  taking 
the  variation  of  the  potential  energy  with  this  assumed  state  and 
setting  the  results  to  zero,  a  condition  of  equilibrium  is 
obtained.  A  load  factor  (eigenvalue)  is  determined  when  the 
second  variation  of  the  potential  energy  ceases  to  be  positive 
definite. 

Imperfection  methods  can  be  used  to  develop  the  equation  of 
equilibrium  for  a  system  with  imperfections.  Eccentricity  of 
load  caused  by  the  imperfection  is  examined,  and  a  solution  is 
obtained.  The  asymptotic  character  of  the  solution  (usually 
radial  displacement)  is  determined.  When  the  deflections  of  the 
structure  increase  without  bound  for  a  relatively  small  increase 
In  applied  load,  a  load  parameter  is  then  said  to  be  determined. 

2.2  Semi-Empirical  Methods 

The  fundamental  approach  or  philosophy  behind  the  employment 
of  semi-empirical  methods  appears  to  be,  in  itself,  quite 
straightforward.  The  classical  buckling  prediction  for  an  ideal¬ 
ized  shell  under  idealized  loads  without  the  influence  of  bound¬ 
ary  conditions  can  generally  be  easily  determined.  For  example, 
the  critical  stress  for  the  axially  loaded  cylinder  can  be 
expressed  as: 

o  =  Z  Et/R  (7) 

cr 

where 

C  =  1.0//3(l-v2) 
or 

C  =  0. 605  for  v  =  0. 3 

The  quantity  C  is  replaced  by  some  functional  quantity  called  a 
knockdown  factor  that  accounts  for  a  proper  fit  of  all  of  the 
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accepted  experimental  data.  For  example,  Weingarten,  et  al.319 
suggests  the  relation: 


C  =  0.606  -  0.546(1  -  e  -  ^  )  (8)  4 

and  shows  that  the  above  equation  fits  the  lower  bound  of  all  the  < 

accepted  data  from  fourteen  different  sources.  The  adequacy  of 
the  semi-empirical  method  is  necessarily  based  upon  a  great  deal 
of  experimental  data  as  well  as  some  correlation  of  this  data  to 
a  specific  formula. 

In  1957,  Gerard  and  Becker^*  112-117  compiled  a  handbook  of 
structural  stability  for  NACA.  Criteria  were  set  for  flat 
plates,  shells,  stiffeners,  etc.,  which  were  based  upon  linear 
energy  methods  supplemented  by  semi-empirical  theory.  Batdorf's 
simplified  methods3®“42  were  used  extensively  in  Gerard's  compar¬ 
ison. 

In  1968,  under  contract  from  NASA,  Baker,  et  al.^1  collected 
various  analysis  procedures  for  shell  structures  encountered  in 
the  aerospace  Industry.  Chapter  3  of  this  shell  analysis  manual 
cataloged  a  variety  of  stability  criteria  that  was  based  upon  the 
concept  of  semi-empirical  methods.  This  work  was  later  condensed 
and  published  by  the  original  authors^2.  Similar  collections 
were  made  by  other  organizations  in  the  aerospace  and  energy 
fields.  Unfortunately,  much  of  this  data  was  considered  proprie¬ 
tary  to  the  supporting  organization. 

With  all  of  the  complex  problems  associated  with  the  deter¬ 
mination  of  a  stability  criteria  for  stiffened  and  unstiffened 
shell  structures,  no  single  criterion  can  be  used.  Bounds  on  the 
application  of  the  criteria  need  to  be  established.  This  was 
done  to  some  extent  by  Baker,  but  it  is  not  complete.  An  attempt 
to  collect  and  computerize  a  semi-empirical  stability  criteria 
was  made  by  Clterley®1.  Criteria  from  Baker^,  Smith  and 
Spier2^8  and  Lakshmikantham  and  Gerard^^  was  used.  The  only 
cons ideratlon  of  imperfection  sensitivity  was  through  a  correla¬ 
tion  factor.  Usually  the  knockdown  factor,  ^  of  Eq.  (8),  would 
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incorporate  this  additional  factor  (y^).  For  example,  the  factor 
y,  given  as: 


Y  =  1 


0. 901 ( 1-e 


-/R7t  /16) 


(9) 


is  suggested  for  curved  isotropic  panels.  Buchert^0  applies  the 
same  basic  idea  to  spherical,  cylindrical  shell-like  structures 
and  reticulated  shell-like  structures  with  positive  and  negative 
Caussian  curvature.  Formulas  are  suggested  that  appear  to  have 
the  same  form  as  the  classical  critical  stress  prediction  but 
considers  an  imperfection  parameter.  A,  that  includes  not  only 
the  imperfect  shape  but  also  the  deformation  of  the  shell  just 
prior  to  buckling.  For  example,  the  critical  buckling  stress  for 
a  spherical  shell  subjected  to  any  loading  is  given  by  Eq.  (7) 
where. 


C  =  -  0.54  ~  -  0.145  H  +  {1.09  (t-)2 
^m 

+  0.359  C^-)3}172 
m 

H  =  [9.9  (£-  )2  +  3.08  C^-)3]172 
m  m 

t  =  effective  membrane  thickness 
m 

t^  =  effective  bending  thickness 

The  cautions  in  adopting  the  knockdown  factor  concept  are  that 
one  is  willing  to  (1)  accept  the  Judgment  of  the  original  author 
as  to  what  is  credible  data  and  (2)  assure  that  the  boundary 
conditions,  loading  and  shell  configurations  are  comparable. 

Miller^5-207  has  collected  a  great  deal  of  information 
using  empirical  data  and  was  the  principal  behind  the  criteria 
suggested  for  Code  Case  N-284  (Metal  Containment  Shell  Buckling 
Methods)  of  the  ASME. 


0.03  H}172  (10) 


19 


At  the  same  time  that  semi-empirical  efforts  were  being 
advanced,  more  experimental  results  became  available.  Instead  of 
establishing  the  precise  behavior  of  shells,  more  questions  were 
being  raised.  The  effects  of  the  testing  machine  compliance, 
probable  boundary  conditions  and,  in  fact,  even  the  buckled 
shape,  were  questioned.  Yoshlmura33^> 335  showed  that  the  cylin¬ 
drical  shell  could  be  developed  from  a  series  of  plane  triangles, 
and  thus  the  shell  could  buckle  into  a  lower  energy  state  (i.e., 
bending  or  folding  along  the  intersection  of  the  triangles).  A 
crude  mathematical  development  of  the  postbuckled  state  corres¬ 
ponds  to  the  assumed  three-term  description  suggested  by  von 
Karman.  Ponsford237  demonstrated  that  the  buckled  pattern  was 
dynamic  and  could  not  be  maintained  even  under  controlled  static 
loads.  Almroth^  found  that  by  increasing  the  number  of  terms 
from  3  to  14,  in  representing  the  radial  displacement  and  the 
initial  imperfections,  the  value  of  the  postbuckllng  load  for  the 
cylinder  would  decrease. 

Further  discussion  of  Donnell's  approach  was  presented  by 
Hoff . 134  whQ  ixiustrated  that  by  evaluating  over  1100  terms  of 
the  total  potential  energy,  that  in  the  limit  the  magnitude  of 
the  radial  displacement  had  a  significant  effect  on  predicting 
the  lower  bound  for  the  postbuckllng  state  and  that  the  Yoshimura 
pattern  could  be  approached.  This  latter  finding  put  a  virtual 
end  to  the  search  of  what  was  then  the  popular  method  of  estab¬ 
lishing  a  lower  bound  in  the  critical  load  through  linear  energy 
methods . 

The  effects  of  boundary  conditions  has  also  been 
reviewed-*-^ » 135, 137, 255, 261, 562, 298#  Ohlra22^  Nachbar  and 

Hoff21^,  Sobel2^,  and  Almro.n*3  also  discussed  the  problem  of 
boundary  conditions.  Almroth  developed  relationships  between  the 
length  of  the  shell  and  edge  restraints  for  eight  possible  bound¬ 
ary  conditions  that  could  exist. 


2.3  Geometric  Nonlinear  Theories 

Continued  research  in  the  adaptation  of  von  Karman's  nonlin¬ 
ear  shell  theory  was  made  by  Leggett  and  Jones and  by 
Michielsen20^.  Unfortunately,  little  information  has  been  pub¬ 
lished  on  the  nonlinear  methods.  Chien^®  presented  one  theory 
that  was  later  criticized  by  Goldenveizer  and  Lur'ye11^  and 
Reiss^^.  The  major  complaint  was  that  Chien's  development  could 
not  handle  displacement  boundary  conditions  correctly. 

pii  1 

Reissner  developed  finite  displacement  equations  of  equilib¬ 
rium  for  the  axlsymmetric  response  of  a  shell  of  revolution. 

This  investigation  was  far-reaching  inasmuch  as  the  extension  of 
linear  methods  had  been  shown. 

The  development  of  the  nonlinear  equations  of  equilibrium 
for  general  shells  has  had  less  intense  study  than  the  linear 
set.  Prom  a  general  point  of  view,  the  discussion  of  Zerna^^ 
gives  the  most  complete  theory  on  the  nonlinear  behavior. 
Naghdi21^  has  also  presented  a  nonlinear  theory.  Perhaps  the 
most  commonly  accepted  theory  is  that  of  Sanders2-^,  which  is 

limited  to  small  strains  and  moderately  small  rotations. 
oho 

Reissner  developed  a  system  of  nonlinear  differential  equa¬ 
tions  for  the  symmetric  deformation  of  shells  permitting  large 
rotations  and  finite  strains.  Danielson^  developed  the  bifur¬ 
cation  and  postbuckllng  equations  using  asymmetric  expansion 
techniques  from  an  energy  criterion.  Yokoo  and  Matsunaga^2 
derived  a  fundamental  set  of  two-dimensional  shell  equations. 

The  structures  under  investigation  were  made  of  rather  simple 
geometric  forms,  although  some  investigation  was  made  of  compos¬ 
ite  vessels,  such  as  the  joining  of  cylinders  and  hemispherical 
heads12^. 


3.0  RECENT  DEVELOPMENTS  IN  ANALYTICAL  METHODS 

As  stated  previously,  the  analytical  methods  used  In  deter¬ 
mining  the  Imperfection  sensitivity  of  shells  usually  follow 
three  avenues: 

(1)  Equilibrium  methods 

(2)  Energy  methods 

(3)  Imperfection  methods 

Each  of  these  methods  has  been  able  to  demonstrate  that  the  crit¬ 
ical  load  of  a  particular  shell  subjected  to  a  given  load  can  be 
reasonably  predicted.  For  many  engineers,  these  predictions  are 
too  cumbersome  and  a  simple  semi-empirical  knockdown  factor  would 
do.  However,  when  the  conditions  under  investigation  go  outside 
established  limits,  the  engineer  should  exercise  caution. 

Further,  the  engineer  must  assess  the  difference  between  the  real 
world  and  the  highly  idealized  buckling  analysis.  In  many 
instances  these  differences  are  collectively  classified  as  imper¬ 
fections.  Some  areas  that  have  been  separately  identified  are: 


(1)  Geometric  imperfections 

(2)  Boundary  conditions 

(3)  Stiffness  or  material  property  variation 

(4)  Nonuniformity  of  load 

(5)  Prebuckling  deformation 

(6)  Dynamic  loads  and  responses 

(7)  Nonlinear  material  behavior 

In  order  to  act  intelligently,  the  engineer  must  be  aware  of 
the  accepted  analytical  methods  and  their  possible  overlapping 
range  of  applicability  and  controversy.  To  do  this,  the  details 
of  each  method  must  be  examined.  Only  the  first  two  methods  will 
be  discussed  in  detail  here.  Because  of  their  complexity,  only 
the  highlights  will  be  presented.  The  derivations  are  exceed¬ 
ingly  difficult  and  generally  require  the  introduction  of  some 


mathematical  hypothesis  that  may  not  be  of  interest  to  the  casual 
reader.  Therefore,  this  presentation  will  only  attempt  to 
demonstrate  the  complexities  associated  with  the  applicable 
method,  and  will  not  provide  a  complete  and  sufficient  derivation 
of  each  method.  The  cylinder  under  axial  compression  has  been 
used  as  the  classic  problem  to  demonstrate  the  conditions  of 
stability  and  the  effects  of  imperfection.  The  discussion  will 
therefore  be  limited  to  this  case. 


3. 1 


Unstiffened  Shells 


3.1.1  Equilibrium  Method 

In  order  to  amplify  the  above  statements,  and  to  gain  a  more 
thorough  understanding  of  the  buckling  branching  paths,  consider 
a  perfect  cylinder  with  a  radius,  R,  thickness,  t,  and  length, 

L.  Assuming  a  radial  displacement,  w,  in  a  shallow  shell 
Cartesian  system  [x,  y  or  Re]  and  with  the  use  of  the  Airy  stress 
function,  $,  the  membrane  stress  resultants  are  defined  by: 


N  = 
x 


yy 


N  =  <b 

y  .xx 


(ID 


N  =  <6 
xy  ,xy 


Donnell's  equilibrium  equations  are  then  defined  in  terms  of  the 
Airy's  function  and  a  displacement  function,  W: 
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Where  V  is  the  biharmonlc  operator  (— ^  and 

X  ^  V 

the  operator  L(Q,V)  is  given  by:  J 


For  a  possible  three  mode  postbuckled  response  denoted  as 

(S  S  and  —  one  axisymmetric  and  two  non-axisymmetric)  with 

L  y  d.  J 

a  cylinder  subjected  to  an  axial  stress,  a,  the  displacement  can 
be  assumed  as: 

w  =  +  W  (14a) 

and  Airy's  function  as: 

4>  =  -  \  aty2  +  $  (14b) 

The  particular  solution  form  of  the  functions  <j>  is  assumed 
to  be  a  double  Fourier  series  such  as  Eq.  (6).  The  more  general 
solution  considers  sine  terms  as  well.  This  particular  solution 
is  substituted  into  Eq.  (12a)  and  the  coefficients  are  determined 
by  equating  like  terms.  Substituting  a  like  function  in  W  also 
including  sine  terms,  if  applicable,  into  Eq.  (12b)  and  applying 
a  Galerkin  procedure,  yields  an  error  function.  This  function  in 
turn  is  used  along  with  orthogonality  conditions  of  the  assumed 
series  to  yield  the  relationship: 

C1(a1-X)  +  (|a2+a3  )U2-?2)  =  0 

f,2[ a4-x  +  8(a2+2a3)  =  0  (15) 

f^Ca^-X  -  8(a2+2a3)  ^  ]  =  0 

Equation  (15)  identifies  the  possible  branching  paths  that  a  per¬ 
fect  shell  can  take  when  multi-modes  exist  as  X  approaches  X* 

(see  Figures  6a  and  6b). 

The  bifurcation  of  a  shell  can  be  reached  under  two 
conditions:  a.^  >  a^  or  a^  <  a^.  For  the  first  case 

(a.^  >  a^),  the  linear  path  I  is  followed  until  x=x#  =  a^.  Then, 
depending  upon  the  participation  of  the  postbuckling  modes,  four 
possible  distinct  paths  can  be  followed: 


^2=S~°  and 
or  ^  *°  and 

^2  or  ^3  and  ^1 

?2  and/or  *0  and  *0 

These  branching  paths  can  be  approached  from  either  the  positive 
or  negative  side  of  each  mode.  When  a^  <  a^ ,  the  paths  are 
slightly  more  restrictive.  First,  path  I  is  followed 
until  X=X*=a^.  Branching  can  then  take  place  when: 


52=C^=0  and  S-j^O  until  5^  = 


a4~al 


then  bifurcation  continues  into  the  mode.  A  similar  response 
in  the  mode  can  take  place  when  <  0.  For  all  of  these 
cases  an  unstable  buckling  path  is  followed. 

The  particular  form  of  the  coefficients  a-L  thru  a^  are  de¬ 
pendent  upon  the  form  of  the  function  w  in  Eq.(l4)  Since  a  sine 
or  cosine  function  is  a  solution  to  Eq.  (12)  w  usually  is  also 
assumed  to  have  the  same  form.  If  an  imperfection  w  is  assumed 
to  exist  with  the  same  content  as  w,  then  the  operators  of 
Eq.  (12)  would  be  4-L(w,w+2w)  and  -L($,w+w)  respectively. 

Arbocz  J  gives  a  rather  complete  description  of  the  above 
procedure  and  Includes  the  effects  of  an  initial  Imperfection  and 
a  three  term  radial  response.  Only  two  simultaneous  modes  of  the 
three  were  considered  in  the  analysis.  Later,  Imbert1-^  examined 
the  three  mode  responses  and  demonstrated  how  the  coefficients  of 
Eq.  (15)  relate  to  the  imperfection  unevenness  parameters 
expressed  by  Donnell  and  Wan.  Both  investigations  amply  demon¬ 
strate  the  mathematical  complexities  they  encounter  when  pursuing 
this  type  of  analysis.  The  fruits  of  their  effort  can  be 
utilized  in  an  efficient  computerized  procedure.  However,  the 
whole  analysis  process  becomes  very  complicated  if  one  tries  to 
include  higher  order  terms  of  the  Fourier  series. 
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3.1.2  Koiter's  Method 

Hotter' s  thesis  on  the  stability  of  elastic  equilibrium  in 
^945169  is  perhaps  the  most  significant  piece  of  work  done  in 
that  decade.  Unfortunately,  his  work  was  not  available  in  Eng¬ 
lish  until  1967^5.  The  original  work  was  performed  between  19^0 

and  1943  but  because  of  World  War  II  his  work  was  isolated  from 
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the  rest  of  the  engineering  community  .  Hotter  formulated  the 

basis  of  the  initial  postbuckling  theory  from  energy  consider¬ 
ed 

ations  which  was  triggered  by  a  hypothesis  suggested  by  Cox^  . 

Hotter' s  statement  of  elastic  stability  assumes  the  exis¬ 
tence  of  a  potential  energy  functional  P[u]  where  u  is  the  dis¬ 
placement  state  which  excludes  inelastic  behavior.  Two  adjacent 
states  are  assumed  to  exist  at  the  bifurcation  load  X*  .  The 
potential  energy  functional  is  expanded  into  a  series  of  like 
terms : 

PX[u]  =  p|[u]  +  PX[u]  +  P*[u]  +  PX[u]  +.  . .  +  eQX[u]  >  0  (16) 

where  the  notation  PX[u]  is  the  collection  of  like  powers  of  mth 
order  of  u  and  its  derivatives  for  a  load,  x,  and  Q^Lu]  includes 
the  effects  of  the  initial  imperfection  and  e  denotes  the  magni¬ 
tude.  By  first  restricting  the  analysis  in  the  neighborhood  of 
the  bifurcation  point,  Kolter  assumed  that  the  displacement  field 
immediately  adjacent  to  the  critical  bifurcation  point  is  essen¬ 
tially  of  the  same  form  as  the  classical  buckling  modes.  Prom 
variational  principles,  it  is  well  known  that  first  variation  of 
the  potential  energy  is  a  necessary  and  sufficient  condition  to 
prove  equilibrium.  Kolter  shows  that  by  expanding  the  potential 
energy  function  for  a  linear  combination  of  adjacent  buckling 
modes,  the  system  is  stable  when  the  second  variation  of  the 
potential  energy  is  semi -definite  positive.  Further  the  third 
variation,  P3[u]=0  and/or  the  fourth  variation  P4[u]>0  for  all 
the  buckling  modes.  The  postbuckling  behavior  of  shells  as 
represented  in  Figure  5  was  suggested  by  Kolter.  Once  the  basic 
idea  had  been  presented,  several  extensions  of  the  method  were 


presented  by  Holter^^-^^  and  others.  Seide2^  presents  a 
rather  condensed,  yet  complete,  description  of  Koiter's  theory 
which  is  applicable  to  shells.  Budiansky  and  Hutchinson  and 
their  colleagues  at  Harvard  have  contributed  to  the  understanding 
of  Koiter's  approach  and  have  successfully  advanced  the  method  in 
several  areas^-^>  1^7-153. 

In  order  to  demonstrate  the  usefulness  of  Koiter's  simple 
yet  elegant  approach,  Budiansky^  in  1967  gave  a  presentation  on 
the  theory  of  initial  postbuckling  behavior  of  shells.  He 
examined  the  asymptotic  behavior  of  the  solution  to  the  governing 
equation  of  a  cylinder  subjected  to  torsion.  Amazigo  and 
Budiansky^  perhaps  have  described  the  more  general  approach, 
while  Stephens2^  applied  the  procedure  to  a  cylindrical  panel, 
and  Fitch10^  applied  the  method  to  spherical  caps.  The  Koiter 
analysis  is  carried  out  by  assuming  that  the  deflection  and 
Airy's  stress  function  can  be  expanded  into  the  form: 


W  =  wo  +  cwx  +  5  w2  +  — 


$  -  $  +  + - 
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(17) 


where  the  subscript  o  refers  to  the  prebuckled  state,  1  and  2  are 
classical  orthogonal  buckling  modes  and  £  represents  the  magni¬ 
tude  of  each  mode.  The  load  factor  X  can  be  expanded  in  terms  of 
powers  of  £  : 

X/X*  =  1  +  a  £  +  b£2  +  -  (18) 


where  x*  is  the  bifurcation  load. 

The  procedure  for  determining  the  coefficients  a  &  b  in  Eq. 
(18)  has  been  outlined  by  Budiansky^.  For  the  symmetric  post- 
buckling  paths  (a=0),  the  parameter  b  for  a  cylinder  can  be 
defined  as: 


b  =  c  //  [2  {*‘,Wj  W^  +  aj'wj  V*2 
-  4*  '(wjw'+W'W^)}  +  *'2' W[)2 
+  $2 *(Wp2  -  24*2^^^]  dxdy  /  // 

where 


a  =  critical  buckling  stress 

c  =  \J 3(  1- v2 i 

Parameter  b  (Koiter's  parameter)  indicates  the  postbuckling  char¬ 
acter  of  the  shell.  When  b  is  negative,  path  II  of  Figure  5b  is 
followed  (unstable).  When  b  is  positive,  path  II  of  Figure  5c  is 
followed  (stable). 

When  a  shape  imperfection  (T  =  w  /t)  Is  introduced,  the 
behavior  of  the  shell  depends  upon  the  amplitude  of  the  Imperfec¬ 
tion.  The  buckling  strength  of  the  axially  compressed  circular 
cylinder  is  greatly  reduced  by  relatively  small  imperfections. 
Kolter  has  shown  that  the  critical  load  factor  A.  is  governed  by: 

(1-A/A*)3/2  =  /-b  (A/A*)  |c|  (20) 

when  a  =  o  and  b  <  o,  as  in  Figure  5b.  The  negative  value  of  the 
factor  b  causes  a  reduction  in  the  buckling  load,  and  is  depicted 
by  the  dotted  line  in  Figure  5b  for  the  imperfect  shell.  A 
restrictive  assumption  implied  In  the  derivation  of  Eq.  (20)  is 
that  £  <  1.  No  attempt  should  be  made  to  apply  this  result  beyond 
this  limit.  Koiter  shows  that  an  axisymmetric  imperfection  can 
cause  a  large  reduction  in  the  load-carrying  capacity  of  a 
cylinder,  as  Indicated  in  Figure  7,  and  he  also  concludes  that 
asymmetric  Imperfections  can  have  an  even  more  pronounced  effect. 
From  the  collection  of  data  from  many  tests,  a  knockdown  factor 


(19) 


dxdy 


29 


between  0.25  and  0.33  has  been  suggested  for  the  axially  loaded 
cylinder.  Kolter  suggests  that  due  to  axisymmetrie  imperfec¬ 
tions,  a  knockdown  factor  of  0.10  can  be  used  as  the  upper  bound 
for  large  imperfections.  However,  Seide  re-examined  Kolter' s 
theory  and  demonstrated  that  the  Kolter  analysis  scheme  is  quite 
limited  and  may  not  accurately  describe  the  failure  mode. 

Most  of  the  investigators  who  have  used  Koiter's  theory  on 

different  shells  of  revolution  have  used  the  buckled  mode  shape 
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as  the  assumed  imperfection  shape.  Lange  and  Newell-10  have 
suggested  that  certain  imperfection  patterns  are  found  on  spheres 
(HEXAGONS)  and  cylinders  (DIAMONDS),  and  Hutchinson-1-^0  has  demon¬ 
strated  that  for  a  spherical  shell  having  a  radius,  R,  and  thick¬ 
ness,  t,  the  Imperfection  shape  could  take  the  form: 

wq  =  Ti  t  cos  (qQ  |)  +  ~l2  t  sin  (qQ  ^)  sin  (qQ  ~  j|)  (21) 

where 


4  - _ , .  2.  /R\2 

q  =  12 ( 1- v  )  (-) 


This  led  to  a  buckling  prediction  equation  of  critical  (p)/clas- 
sical  (p*)  pressure  ratios: 


(1  -  p/p*)2  -  U2I  P/P*  ;  =  °,  ?2  *  0 

(1  -  p/p*)2  =  9 ^1~v-  }  ?!  p/p*  ;  ?!  >  0,  I2  -  ° 


These  same  equations  were  verified  by  Reissner2^  using  certain 
orthogonality  relations  in  Koiter's  procedure  "as  a  consequence 

Op 

of  the  nonappearance  of  secular  terms."  Clterley,  et  al.  ,  has 
demonstrated  through  Imperfection  methods  that  when  the  mode 
shape  is  comprised  of  an  axisymmetrie  and  non-axisymmetric  imper¬ 
fections  (^!  >  0;  C2  >  0)  a  further  reduction  is  realized.  More¬ 
over,  the  influence  of  a  boundary  interacts  with  any  shape  imper¬ 
fection  to  reduce  the  buckling  load. 


Since  shell  structures  are  fabricated,  they  are  subject  to 
manufacturing  tolerances.  In  an  attempt  to  predict  the  buckling 
behavior  of  these  shells,  it  must  be  assumed  that  the  imperfec¬ 
tions  occur  in  a  random  manner^.  For  cylinders,  Fersht101*  found 
that  axlsymmetric  random  imperfections  reduced  the  buckling  load 
less  than  nonaxisymmetr ic  random  imperfections.  Other  imperfec¬ 
tion  geometries  have  been  studied  to  determine  what  local  effects 
of  assumed  shapes  have  on  the  overall  stability  characteristics 
of  a  shell.  Some  of  the  most  interesting  comparisons  with  exper¬ 
imental  results  in  this  area  are  those  presented  by  Tennyson  and 
his  associate  Muggeridge^^ ,  Caswell, Chan^91~296  and 
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Hutchinson  J  . 

?? 

Finally,  Arbocz,  in  his  thesis  ,  demonstrated  than  an 
adequate  prediction  of  the  knockdown  factors  can  be  performed  if 
the  shape  of  the  imperfection  is  accurately  known.  In  subsequent 
analyses  the  effects  of  boundary  conditions  were  studied  dome 
attempts  to  store  the  imperfection  data  of  test  specimens  into  a 

pO 

data  bank  have  also  been  implemented  . 

3.1.3  Thompson's  Extension  of  Koiter's  Method 

One  group  that  has  developed  a  comprehensive  understanding 
of  Koiter's  principles  is  at  University  College,  London,  and 

consists  of  Thompson  and  Hunt^99-302^  Roorda^^  and  Supple*^?- 

qOQ 

r-°°.  It  has  been  demonstrated  by  Roorda  that  certain  structures 
could  behave  in  a  manner  that  is  not  necessarily  predicted  either 
by  the  equilibrium  or  the  imperfection  methods.  Thompson  and 
Hunt-^^  have  given  a  rather  complete  mathematical  description  of 
Implementing  conceptually  the  energy  theorems  to  an  optimum 
i t  rue tu re . 

Supple”^  illustates  the  basic  concept  of  Thompson's 
generalized  coordinates  in  the  general  theory  of  elastic 
stability  through  an  example  of  symmetric  structural  systems. 

The  potential  energy,  P,  of  a  structural  system  is  assumed  to  be 


•  i  Curie  I,  Lon  of  a  genoral  Load  parameter  1'  and  two  non-dlinen- 

s Lonal  generalized  coordinates,  and  .  The  latter  two 
represent  possible  buckling  modes.  For  an  incremental  change  in 
load,  6f ,  from  an  inital  value,  fQ,  the  potential  energy  can  be 
written  in  the  form: 


P(f,ui)  =  P(fQ  +  fif^) 


Expanding  the  right  hand  side  in  the  form  of  a  Taylor  series 


Ptf.up  -  P(fo,0)  +  §  u,  +  «  if 

p  1  o  o 

2  p 
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+  JT  3u.  3u.  Su^  r  uiu  JUk  +  3  3u .  3u,3f  „  UiUj6f 
L  j  o  1  J  fo 


where  the  subscripts  i,j,k  are  able  to  take  the  value  of  1  and  2. 

If  the  buckling  modes  u-^  and  u2  are  assumed  to  be  small,  the 
two  equilibrium  equations  for  a  symmetric  system  take  the  form: 
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Three  solutions  exist  for  a  doubly  symmetric  system.  Two 
solutions  referred  to  as  uncoupled  modal  solutions  are  given  by: 


33P 

3  !  (  5f  -I  )  -V- 
3u^3f 


Up  =  0 


(26) 
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u^  =  0 


(27) 


where  6f^  and  Sf2  are  different  critical  load  values;  i.e.,  at 
ul=u2=0*  fl  ’  fo  +  5fl  •  f2  ’  fo  +  5f2 


corresponding  to  the  different  buckling  modes  of  u^  and  u2.  The 

quadratic  character  of  these  solutions  follows  the  individual 
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bifurcation  paths  in  Figures  5b  and  8a  for  3  P/3u.  <  0,  and 

o  4  4  1 

Figures  5c  and  ob  for  3  P/au^  >  0.  Since  i  represents  two 

possible  modes,  two  additional  combinations  can  exist: 


4  4 

a  p/au^  <  o 


4  4  4  4 

a  p/au2  >  0  and  a  P/auJ  >  0 


4  4 

34P/3u2  <  0, 


These  are  illustrated  in  Figure  8c. 

The  third  solution  is  for  the  coupled  mode,  which  is  satis¬ 
fied  by: 


au?  af 
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(28a) 


where  Af  =  f^  -  f2 

Four  possible  bifurcation  conditions  can  exist  for  the  coupled 
modes : 

Condition  1:  A  transition  path  between  the  uncoupled  post- 
buckling  equilibrium  paths. 

Condition  2:  Always  from  the  primary  uncoupled  postbuckling 
path. 


t  it! 


Condition  3:  Always  from  the  secondary  uncoupled  postbuck- 
ling  path. 

Condition  4:  Prom  one  or  both  of  the  uncoupled  modes. 


When  both  uncoupled  modes  are  rising  (Figure  8a),  the  coupled 
buckling  path  from  the  primary  equilibrium  paths  will  be  rising 
(symmetric-stable)  when: 


34P  > 


(28b) 


and  will  be  falling  (symmetric-unstable)  when: 


2  I 
3u1  3u2 


(28c) 


When  either  or  both  uncoupled  postbuckling  paths  are  unstable, 
the  coupled  buckling  paths  branching  from  the  primary  uncoupled 
path  will  always  be  unstable.  When  one  or  both  uncoupled  paths 
are  stable,  the  coupled  path  from  the  secondary  path  will  always 
be  stable.  When  both  modes  are  unstable  (Figure  8b),  branching 
from  the  secondary  path  will  be  stable  when  Eq.  (28b)  is  met,  and 
unstable  when  Eq.  (28c)  is  met.  The  complexities  of  evaluating 
the  stability  characteristics  of  a  shell  when  employing  this 
extension  of  the  Koiter  method  begins  to  overshadow  the  final 
results . 

In  Koiter' s  method,  a  single  parameter,  b,  defined  by  Eq. 
(19),  is  used  for  the  determination  of  a  stable  shell.  In 
Thompson's  approach,  rather  than  developing  a  single  parameter,  a 
more  general  approach  is  used.  Through  the  employment  of  calcu¬ 
lus  of  variation  principles,  different  conditions  of  the  total 
strain  energy  can  be  examined.  It  appears  that  Thompson's  exten¬ 
sion  of  Koiter' s  method  not  only  gives  the  same  results  as  the 
original,  but  also  predicts  some  additional  behavior.  A  compar¬ 
ison  of  Fig.  6a  and  Figure  9lii  suggests  the  similarity  of  the 
two  methods. 
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Figure  9.  Forms  of  coupled  postbuckling  for  ideal 
doubly-symmetric  structural  systems 
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Figure  ]  0 .  Forms  of  postbuckling  equilibrium  path  for 
doubly-symmetric  structural  systems  with 

imperfections  e  ( -  stable  equilibrium  path 

-  -  unstable  equilibrium  path)  [287]. 


The  important  finding  from  Thompson’s  approach  is  that  even 
with  two  uncoupled  buckling  modes,  which  are  on  stable  rising 
paths,  unstable  coupled  paths  are  possible.  This  phenomena  would 
have  been  ignored  by  other  methods.  Similar  conclusions  are 
drawn  for  systems  with  imperfections  (Figure  10). 

In  a  more  abstract  approach,  Thompson^^  has  utilized 
catastrophe  theory  in  order  to  establish  the  precise  number  and 
type  of  imperfections  to  which  a  system  will  be  most  sensitive. 

He  suggests  that  instead  of  a  two  mode  interaction  representation 
of  imperfections,  seven  control  parameters  may  be  required;  one 
of  which  is  the  load.  Thus,  six  control  parameters  are  then 
required.  The  surface  that  represents  one  of  the  fundamental 
control  parameters,  a  hyberbolic  umbilic,  as  described  by 
ilunt^1,  is  illustrated  in  Figure  11.  Although  this  approach  has 
not  been  thoroughly  tested  or  accepted,  it  may  suggest  why  propo¬ 
nents  of  one  or  another  type  of  imperfection  parameter  may  be 
partially  correct  in  their  findings. 

Although  the  details  may  be  too  idealized  for  many,  two 
quite  simple  theorems  are  postulated  by  Thompson  and  Hunt  that 
tend  to  lay  a  foundation  for  the  understanding  of  bifurcation  in 
shell  structures  and  the  mechanisms  that  induce  the  failure: 

"Theorem  1:  An  initially  stable  (primary)  equilibrium  path 
rising  monotonically  with  the  loading  parameter  cannot 
become  unstable  without  intersecting  a  further  distinct 
(secondary)  equilibrium  path. 

"Theorem  2:  An  initially  stable  equilibrium  path  rising 
with  the  loading  parameter  cannot  approach  an  unstable 
equilibrium  state  from  which  the  system  would  exhibit  a 
finite  dynamic  snap  without  the  approach  of  an  equilibrium 
path  (which  may  or  may  not  be  an  extension  of  the  original 
path)  at  values  of  the  loading  parameter  less  than  that  of 
the  unstable  state." 
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These  theorems  would  appear  to  be  a  rather  complicated 
statement  for  such  a  simple  mechanism  of  single  path  bifurcation. 
On  the  other  hand,  for  the  more  general  case,  complex  functionals 
may  be  required  to  represent  the  total  energy.  This  means  that 
more  complicated  variational  methods  need  be  applied  to  make  a 
proper  evaluation.  However,  with  these  theorems  in  mind,  a 
clearer  understanding  is  developed  on  the  possible  paths  or 
branching  mechanisms  of  shell  buckling  into  symmetric  and  asym¬ 
metric,  stable  and  unstable  postbuckling  modes.  This  type  of 
evaluation  can  best  be  done  using  energy  methods.  Since  it  is  a 
nonlinear  problem,  the  solution  to  a  specific  set  of  conditions 
is  not  unique  and  is  path-dependent.  The  path  that  the  solution 
follows  has  been  shown  to  be  influenced  by  Imperfections.  Energy 
methods  can  examine  the  overall  behavior  of  a  shell's  buckling 
characteristics  and  establish  the  final  equilibrium  position  (if 
one  exists). 


3. 2  Imperfection  Sensitivity  of  Stiffened  Shells 

A  great  deal  of  effort  has  been  expended  in  recent  years  in 

an  attempt  to  more  fully  understand  the  instability  of  stiffened 
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shells.  Examples  are  Amazigo  >  ,  Brodsky-^’^  ,  Bronowicki"^  , 

lard75,  Goldberg118,  Heard130,  Khot163’164,  Lee190,  McGinley201, 
Pederson4’'"9,  Singer271,  Tvergaard307  and  Viswanathan312 .  Because 
of  the  introduction  of  semi -rnonocoque  design  in  aircraft  struc¬ 
tures  and  the  recent  employment  of  integral  stiffener  designs  in 
rocket  and  missile  structures,  interest  has  shifted  to  more 
economical  and  efficient  stiffened  shell  structures.  Again,  when 
optimisation  procedures  for  weight  reduction  are  considered,  the 
stability  of  the  shell  structure  becomes  important.  As  with 
u ns t i p fened  shell  research,  the  effects  of  imperfections  on  the 
performance  of  the  shell  structure  became  one  of  the  key  para¬ 
meters  to  study. 
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3.2.1  Equilibrium  Methods 

Prior  to  the  1960’s  the  popular  approach  in  stability  analy¬ 
sis  of  a  stiffened  shell  structure  was  to  examine  the  response  of 
the  shell  as  two  components^*  289,308,309.  The  stringer  and  ring 
components  that  are  comprised  of  beam  elements  were  analyzed 
separately  from  the  skin.  It  was  assumed  that  the  stringers 
would  carry  the  axial  thrust  and  bending,  while  the  skin  would 
transmit  only  inplane  shear.  Perhaps  the  shear  flow  method  of 
analysis  had  something  to  do  with  this  assumption.  Later,  in  an 
attempt  to  more  accurately  define  the  response  of  the  two  com¬ 
ponents,  the  membrane  capacity  of  the  shell  was  included  in  the 
response,  and  the  overall  response  of  the  shell  was  examined. 

As  the  number  of  stiffening  elements  increases,  one  popular 
method  that  has  been  used  is  to  alter  the  properties  of  the  shell 
by  "smearing"  the  stiffener  properties  into  a  combined  ortho¬ 
tropic  shell.  This  procedure  has  some  drawbacks  because  one  may 
lose  sight  of  the  "thinness”  assumption  of  the  shell  Inasmuch  as 
the  effective  thickness  would  Increase,  which  in  turn  would  tend 
to  increase  the  effects  of  transverse  shear.  Further,  when  con¬ 
sidering  the  stability  of  the  reinforced  shell,  the  discreteness 
of  the  stiffeners  may  tend  to  play  an  important  role.  With  this 
approach  fully  explored,  the  next  extension  would  be  the  fully 
coupled  problem. 

In  19^7,  van  der  Neut^10  determined  that  for  longitudinally 
and  circumferentially  stiffened  cylindrical  shells,  the  eccentri¬ 
city  of  the  stiffener  had  primary  importance  —  inside  stiffening 
reduced  the  buckling  load  up  to  one  third  of  the  case  without 
eccentricity.  (The  greater  efficiency  of  outside  stiffening, 
though  apparent  from  the  general  formula,  was  not  explicitly 
mentioned  by  van  der  Neut. )  This  work  was  done  in  19^2  and  was 
meant  to  apply  to  aircraft  fuselages. 

Later,  Singer,  et  al.^>265  presented  a  similar  finding  — 
that  outside  stiffeners  were  more  efficient  than  inside  stif¬ 
feners,  and  centrally  placed  stringers  were  the  least  efficient. 
Nash  7 ,  using  small  deflection  theory,  examined  ring  stiffened 


shells  under  external  pressure  such  as  used  in  submarines. 
Bijlaard^  accounted  for  the  torsional  restraint  of  the  rings 
under  the  same  loading  condition.  Becker^  used  the  linearized 
Donnell  equations  for  the  stability  analysis  of  stiffened  cylin¬ 
ders  subject  to  axial  compression,  torsion  and  hydrostatic  pres¬ 
sure  . 

In  I960,  Thielemann^?  developed  a  nonlinear  theory  for  the 
buckling  of  orthotropic  shells.  He  examined  the  nonlinear  post- 
buckling  behavior  of  cylinders  and  considered  an  infinitesimal 
imperfection  and  radial  displacements,  and  ignored  boundary  con¬ 
ditions.  Through  a  Galerkln  procedure,  the  classical  buckling 
load  was  defined  for  cylinders  with  various  stiffening  character¬ 
istics.  By  examining  the  effect  of  these  stiffening  characteri¬ 
zations,  Thielemann  determined  that  the  resulting  deformation 
patterns,  especially  those  that  had  short  wave  lengths  in  the 
longitudinal  and  ci rcumf erential  directions,  thus  forming  a 
diamond  shape  pattern,  was  strongly  Influenced  by  the  direction 
of  stiffening. 

The  key  parameter  that  Thielemann  had  used  in  his  assessment 
is  r,  which  is  the  ratio  of  the  products  of  bending  and  membrane 
stiffness  in  each  direction.  For  an  isotropic  unstiffened  shell, 

Y  -  1.  If  a  cylinder  is  stiffened  by  longitudinal  members, 

Y  '  1 •  For  ring  stiffened  cylinders,  y  <  1.  The  typical 
behavior  of  an  axial  stiffened  and  circumferential  stiffened 
she'll  subjected  to  an  axially  compression  is  depicted  in  Figures 
1?  and  13.  As  seen  in  a  comparison  in  Figure  12a  and  Figure  12b, 
a  dramatic  reduction  in  the  load  parameter  X  (normalized  to  the 
classical  value  X*)  is  experienced.  Similar  reductions  are 
experienced  in  the  end  shortening,  6  and  the  classical 

value  6  .  The  effects  of  axial  stiffening  are  similar  to  that 

of  internal  pressure,  while  for  circumferentially  stiffened 
shells,  only  a  slight  effect  is  experienced  (Figures  lA  and  15). 

Thielemann  has  shown  by  normalizing  with  extensional  (Apj) 
and  bending  (D|j)  stiffnesses  that  the  classical  buckling 
load,  N  ,  connected  with  a  checkerboard  buckling  pattern  is 


" igure  15.  Postbuckling  functions  of  axially  compressed  pressurized 
orthotropic  cylinders  -  circumferentially  stiffened. 


increased  appreciably  when  the  internal  pressure  P  is  increased. 

The  ring  buckling  load,  N  ,  is  the  upper  limit.  At  the  points 

c*o 

characterizing  the  classical  loads,  the  postbuckling  functions 
corresponding  to  diamond-shape  buckling  bifurcate,  and  for  large 
values  of  the  end  shortening  e  ,  they  asymptotically  approach  the 
ring  buckling  load. 

Smith,  Spier  and  Fossum2"^  generated  a  number  of  parametric 

studies  for  an  unpressurized  stiffened  cylindrical  shell  with 

clamped  boundary  conditions  based  upon  Thielemann's  method.  A 
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similar  study  was  made  by  Milligan,  et  al.  ,  for  axial 
compression,  torsion  and  hydrostatic  pressure.  Almroth^  assumed 
a  five-term  displacement  function  for  the  radial  displacement  and 
re-evaluated  Thielemann's  method  for  different  extensional  and 
bending  rigidities.  Later,  the  influence  of  stiffener 
eccentricities  were  examined  through  the  coupling  of  the  membrane 
and  bending  and  twisting  coefficients  of  orthotropic  shells  °. 

Discreteness  of  ring  stiffeners  was  investigated  by  van  der 

a  i  n 

NeutJ  for  the  cases  of  short  longitudinal  wave  length.  He 
found  that  the  error  in  the  buckling  load  for  smeared  rings  is 
only  about  one  percent  too  great  when  the  average  half  wave 
length  Is  greater  than  twice  the  ring  pitch.  With  smaller  wave 
lengths,  the  error  increases  —  a  five  percent  error  of  a  one  and 
one  half  ring  pitch.  Block^-^2  later  investigated  the  effects 
of  eccentricity  and  discreteness  of  an  eccentrically  stiffened 
shell  and  found  that  the  predicted  buckling  load  may  be 
substantially  affected  by  considerations  of  prebuckling 
deformations ,  eccentric  loading,  and  discreteness  of  ring 
stiffeners . 

One  of  the  major  efforts  in  understanding  the  response  of 
stiffened  shells  has  been  performed  at  the  Israel  Institute  by 
Singer2^-2?5  and  associates,  Baruch-^-^ ,  Haftka-1-2^, 
Rosen2^~253  aruj  Weller^2^-^2^ .  Over  thirty  reports  have  been 
generated  on  the  theoretical  and  experimental  verification  of 
buckling  of  cylindrical  shells.  Stringer  construction  with  non- 


uniform  cross-sections  were  examined.  A  rigorous  analysis  that 
takes  into  account  the  discreteness  of  the  stiffener  was  per¬ 
formed.  Unlike  Block’s  findings.  Singer  shows  that  unless  the 
number  of  rings  is  less  than  three  or  four  for  a  shell  with  L/R  > 
2,  the  difference  between  the  "smeared"  and  "discrete"  theories 
is  negligible. 

The  analysis  procedure  used  by  Singer,  et  al.,  is  to  examine 
the  variations  of  the  equations  of  equilibrium  in  terms  of  dis¬ 
placements,  (i.e.,  equilibrium  methods).  The  stress  and  moment 
resultants  are  given  in  terms  of  the  displacements  for  the  cir¬ 
cumferential  and  longitudinal  harmonic. 

Through  a  Galerkin  method,  a  system  of  determinates  results 
in  a  tridiagonal  symmetric  eigenvalue  problem.  Singer  claims 
that  the  accuracy  of  the  solution  to  these  problems  is  quite 
good,  even  up  to  a  matrix  order  of  100;  however,  the  effects  of 
an  imperfection  is  not  examined. 

3.2.2  Koiter’s  Energy  Methods 

Brush^-59  represents  the  stability  problem  of  an  ortho- 
tropically  stiffened  cylindrical  shell  as  the  statement  of  two 
adjacent  equilibrium  states  —  one  unbuckled,  and  the  other  in  a 
buckled  state.  The  solution  procedure  follows  Koiter's  approach 
by  stating  the  potential  energy  increment  in  terms  of  the  Euler 
equations  for  the  second  variation  of  the  potential  energy.  The 
restriction  is  that  only  simply  supported  shells  are  considered, 
although  the  torsional  rigidity  of  the  stiffeners  are  included. 
The  solution  is  obtained  using  a  finite  difference  representation 
of  the  energy  increment.  Stephens^^  examined  the  imperfection 
sensitivity  of  cylindrical  panels  and  accounted  for  the  torsional 
rigidity  of  a  longitudinal  stringer.  Again,  Koiter’s  aproach  was 
used. 

in  a  series  of  studies,  Budiansky^’  Hutchinson1^1  >  *52 
and  .Stephens'1^  examined  the  postbuckling  behavior  of  stiffened 
cylindrical  shells.  The  effects  of  eccentricities,  prebuckling 
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deformations,  imperfections,  boundary  conditions,  pressure  and 
other  loading  conditions  were  considered.  In  all  studies  pre¬ 
sented  by  this  group,  the  orthotropic  properties  were  based  upon 
a  "smeared”  representation  of  the  stiffeners. 

Hutchinson  and  Frauenthal^3  developed  a  postbuckling  analy¬ 
sis  for  orthotropic  shells  that  parallels  that  presented  for  a 
uniform  shell  represented  by  Eqs.  (12)  through  (19).  The 
development  involved  the  Donnell  equations  written  in  the  form: 


For  compatibility  - 


L„[F]  -  Ln  [W ]  -  W  +  W  W  +  W  W 
H  Q  , xy  ,xx  , yy  ,yy  ,xx 


(29) 


-  2W,xy  W,xy  -  0 


For  equilibrium  - 


LrlW]  +  L~[F]  -  FW  -FW  +2FW 
D  QL  J  ,xx  ,yy  ,yy  ,xx  ,xy  ,xy 


(30) 


-F  W  -F  W  +2F  W  +  p  =  0 
,xx  ,yy  ,yy  ,xx  ,xy  ,xy 


where  the  operators  LH,  Lq,  Ld  are  functions  of  the  shell, 
stringer  and  ring  properties. 

The  above  equations  are  in  terms  of  the  Airy  stress  func¬ 
tion,  F,  the  radial  displacement,  W,  the  imperfection 

A 

func.tLon,  W,  internal  pressure,  p,  shell  properties  (thickness, 
modulus  of  elasticity  and  Poisson's  ratio)  and  stiffener  proper¬ 
ties.  For  a  given  set  of  boundary  conditions  and  a  state  of 
stress,  three  solutions  of  the  W  -  F  function  must  be  determined. 


It  is  assumed  that  each  function  has  the  form: 


1  r  voR  ,  *  . 

W  =  — g—  +  w  +  w 


(3D 


1  2  * 

F  =  -  5-  ohy*  +  F  +  f 
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The  first  term  Ln  each  function  corresponds  to  the  membrane 
state.  The  starred  terms  correspond  to  the  prebuckled 
state  (W^°\  F^);  and  the  last  term  corresponds  to  the  Initial 
postbuckling  stage  and  is  comprised  of  two  orthogonal  states  - 
W(1),  W(2)  and  F^1*,  F  (2). 

Again  following  Koiter's  basic  approach,  each  function  can 
be  expanded  about  the  classical  load  in  the  form  of  Eqs.  17  and 
18.  For  a  symmetric  system,  in  Eq.  18, 

a  =  0 


F(2)*(W(1)aW(1) 


+  2F(1)HW(1)  Jf(2) 


P  [F°*(W(1),W(1) )  +  2F(1)*(W°,W(1) )] 
o  c  c 

wLth  the  notation: 

A* ( B , C )  =  /  [A  BC+A  B  C 

Js  ,xx  ,y  ,y  ,yy  ,x  ,x 

-  A  (B  C  +  B  C  )]ds 
,xy  ,x  ,y  ,y  ,x'J 


•  >  -  8(  ) 
>c  ~  a  p 


It  follows  that  an  imperfection  of  amplitude  6  will  define  the 
load  factor  just  prior  to  buckling.  Koiter  has  shown  that  the 
value  of  b  determines  the  stability  of  the  shell,  per  Eq.  (20) 
for  cylinders,  repeated  here  as: 

X  3/2  1/2  -  X 

(1  -  -jS-)  -  (-B)  |i|  _£  (3 

c  c 

with  b  replaced  by  b. 

When  the  prebuckling  state  is  pure  membrane  b  =  b,  but  when  pre 
buckling  nonlinear  states  may  be  important: 


—  u  I 


/,.,(!)  o  „(1K  1 2 


^ 1 ) ,  W  ^ 1 '  )  +  2F 


o  ,,  1) 


After  much  algebra^**,  the  functions  W  and  F  were  evaluated 
for  a  simply  supported  shell.  For  the  axially  stiffened  cylinder 
under  axial  compression,  Hutchinson  illustrates  comparative  stif¬ 
fening  effects.  The  outside  stiffened  cylinder  denoted  by  its 

eccentricity  e  is,  in  general,  more  imperfection-sensitive  than 
s 

one  with  inside  stiffeners.  The  buckling  load  may  be  well  below 
the  classical  value  for  light  and  heavy  stiffeners  (denoted  by 
the  extensional  ratios  of  cross-sectional  area  As  to  panel  width 
dg  and  stiffener  bending  rigidity  EIS  to  plot  flexural  rigidity, 
dgD);  see  Figures  16-18.  Also,  the  effect  of  stringer  eccentri¬ 
city  is  less  prominent  under  hydrostatic  pressure  than  under 
axial  compression  (same  as  Thielemann's,  van  der  Neut's  and 
Singer's  conclusions).  The  same  is  true  for  ring-stiffened 
shells.  Hutchinson  also  concludes  that  "if  one  takes  into 
account  the  predicted  insensitivity  of  the  axially  stiffened 
cylinder  in  the  lower  range  of  Z  (Batdorf's  parameter)  and  the 
sensitivity  of  the  ring-stiffened  specimens,  then  the  advantage 
of  axial  stiffening  is  even  more  pronounced.” 

Several  studies  have  been  made  to  examine  the  imperfection 
sensitivity  of  stiffened  shells  when  overall  and  local  buckling 
modes  couple.  Byskov  and  Hutchinson^  suggest  an  approach 
extending  Kolter's  method.  Koiter1^  demonstrated  that  Byskov 
and  Hutchinson,  in  a  preliminary  report,  ignored  a  crucial  term 
in  the  mode  interaction  which  had  an  important  effect  on  the 
result.  This  criticism  brings  out  one  of  the  major  disadvantages 
of  the  asymptotic  methods  —  namely,  the  importance  of  the  terms 
omitted  in  the  expansion  cannot  always  be  measured  unless  the 
system  is  resolved  with  the  greater  complexity  included. 

Koiter's  method  calculates  the  slope  and  curvature  of  the  post- 
buckling  path  at,  or  immediately  adjacent  to,  the  bifurcation 
point.  In  some  cases,  the  asymptotic  representation  of  the  post- 
buckling  path,  based  upon  these  two  quantities,  may  be 
misleading;  hence,  imperfection  sensitivity  results  might  be 
asymptotically  valid  but  practically  misleading. 
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The  Imperfection  sensitivity  of  stiffened  shells  is  farther 
complicated  by  the  fact  that  multi-mode  interaction  tends  to  mask 
the  effects  of  the  individual  components  on  individual  or  single 
mode  buckling  behavior.  Tvergaard^^  reviews  this  area  and  sug¬ 
gests  that  single  mode  asymptotic  techniques  are  inaccurate  for 
imperfection  sensitivity  predictions.  van  der  Neut^H  shows  that 
for  built  up  columns  where  the  separate  local  and  overall  buck¬ 
ling  occurs  at  the  same  stress,  the  interaction  of  the  two  modes 
causes  unstable  equilibrium  with  the  effect  that  local  imperfec¬ 
tions  yield  a  reduction  of  the  overall  buckling  load.  This 
result  was  confirmed  by  Koiter  and  Kuiken^^.  Koiter  and 
Pignataro^S  examined  stiffened  panels  and  concluded  that  only  a 
small  imperfection  sensitivity  should  be  realized.  Koiter  and 
van  der  Neut1^  show  that  the  earlier  analysis1^  is  Inadequate 
for  panels  where  in  the  local  mode  stiffener  wall  deflections  and 
skin  deflections  have  an  equal  order  in  magnitude.  The  first  two 
local  modes,  together  with  the  overall  mode,  are  needed  for 
assessing  the  effects  of  the  Interaction. 


4.0  NUMERICAL  METHODS 


There  have  been  some  attempts  to  evaluate  the  imperfection 
sensitivity  of  shells  numerically,  particularly  for  shells  of 
revolution.  Cohen^3“86  used  Koiter's  unique  mode  bifurcation 
buckling  theory  for  the  axisymmetric  nonlinear  prestressed  states 
of  a  ring-stiffened  shell  of  revolution  with  an  arbitrary  imper¬ 
fection  shape.  Haftka,  et  al.^^,^3,  used  the  modified  struc¬ 
ture  method  for  two-dimensional  finite  elements  and  found  some 
limited  success  in  this  approach.  They  found  that  the  nonlinear 
states  of  stress  prior  to  buckling  caused  buckling  to  occur  at  a 
load  significantly  different  from  that  obtained  through  linear 
stability  predictions.  Maewal  and  Nachbar1^  examined  the  post- 
buckling  behavior  of  a  clamped  uniform  cylindrical  shell  under 
uniform  axial  compression  using  the  same  method.  They  required  a 
high  order  isoparametric  finite  element  to  gain  the  accuracy 
required  for  imperfection  sensitivity  predictions. 

The  existence  of  closely  spaced  bifurcation  points  leads  to 
a  complex  and  sometimes  time-consuming  search  for  the  minimum 
eigenvalue.  Again,  the  theories  of  Thompson  seem  to  apply.  If 
the  analysis  procedure  would  permit  response  in  all  Fourier  har¬ 
monics,  then  a  minimum  imperfection-sensitive  calculation  could 
be  acornpllshed.  If  discrete  response  Fourier  harmonics  or 
limited  models  are  used,  then  the  minimum  response  harmonics  and 
all.  important  coupled  modes  must  be  known  a  priori. 

The  use  of  the  direct  numerical  evaluation  of  the  nonlinear 
behavior  of  shells  has  been  attempted  with  some  limited  success. 
Several  computer  programs  exist  that  are  capable  of  examining  the 
static  imperfection  sensitivity  of  shells  —  BOSOR^^,  STAGS?, 
FASOK^?,  and  NBALL (S ATANS ) .  The  last  three  permit  general 
imperfections ,  while  BOSOR  permits  only  axisymmetric 
imperfection.  Only  STAGS  and  NBALL  permit  general  loading  to  be 
considered.  STAGS  can  examine  the  general  two-dimensional  shell 
form,  while  NRALL  is  limited  to  shells  of  revolution. 
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There  have  been  some  special  programs  that  have  been  written 
In  an  attempt  to  reduce  the  computer  run  time  for  a  geometric 

ho 

nonlinear  analysis.  Bauld  and  Satyamurthy  J  have  written  a  pro¬ 
gram  based  upon  the  finite  difference  energy  method  (similar  to 
that  used  in  STAGS)  for  a  cylindrical  panel.  The  examination  of 
the  sign  of  the  final  determinant  is  used  as  a  key  for  the  limit 
(bifurcation)  point.  Since  a  Newton-Raphson  procedure  is  used, 
only  the  static  nonlinear  prebuckled  states  up  to  buckling  can  be 
studied.  To  examine  beyond  prebuckled  states,  the  dynamic  option 

must  be  used.  Only  STAGS  and  NBALL  have  this  capability. 

1  6S 

Kicker  J  has  written  a  computer  program  that  examines  the  linear 
bifurcation  predictions  for  a  cylinder  with  longitudinal  and 
circumferential  stiffeners  or  panels  with  edge  stiffeners  using  a 
selected  theory  of  Plugge,  Timoshenko  or  Donnell.  Smeared  stif¬ 
fener  properties  are  assumed,  and  eight  possible  boundary  condi¬ 
tions  are  available. 

An  interesting  comparison  between  Koiter's  initial  postbuck- 
llng  theory  and  numerical  results  from  the  nonlinear  analysis 
program  NBALL  has  been  made  by  Stillwell  and  Ball2^.  From  Fitch 
and  Budiansky10^  the  initial  slope  of  the  bifurcation  branch  of 
the  equilibrium  path  on  a  load-displacement  plot  had  been  pre¬ 
dicted  for  a  uniformly  loaded  spherical  cap  with  a  rise  height 
parameter  of  A=8.  The  numerical  results  were  obtained  by 
assuming  a  slight  eccentricity  of  load  of  the  form: 

p  =  Pq  (  1  +  e  cos  n  0  ) 

3y  first  obtaining  the  minimum  critical  Fourier  harmonic,  n, 
Stillwell  and  Ball  determined  the  maximum  response  for  the  shell 
under  different  increasing  eccentricities  to  coincide  with 
Fitch’s  predictions.  Since  the  eccentricities  were  small,  two 
important  conclusions  can  be  drawn: 

l)  Numerical  methods  can  be  used  to  accurately  predict  the 
bifurcation  loads  of  shells  of  revolution  when 

imperfection  sensitivity  is  of  concern. 

Koiter's  method  is  applicable  for  small  imperfections. 
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5.0  EXPERIMENTS 

The  ultimate  test  of  a  particular  theory  or  numerical  method 
is  the  comparison  of  predictions  to  experimental  data.  Without 
the  full  understanding  of  the  static  stability  behavior  of 
shells,  experiments  were  performed  with  limited  success.  Prom 
the  early  tests  of  Robertson2^ ,  Donnell^,  Tokugawa^®^ , 
Lundquist'*' and  others  such  as  Bridget-’**,  Holmquist^^, 
Ikeda1'’^,  Kanemitsu1^0,  and  at  NACA210, 226, 282  in  1930*3,  and 

by  NACA-sponsored  tests  in  the  1940  * s 12-18 , 77 , 92 , 93 , 101 , 102 , 132 , 
140,141  inconsistent  data  was  reported.  Since  the  cylinder 
possesses  the  greatest  sensitivity  to  axial  load,  this  shape  was 
used  for  most  of  the  experimental  studies  in  the  50*  s  and  60  *  s  at 
NASA23.M.7'M00,209,231,23‘l,2M,285  and  5y  oalletly110. 111 , 

Allen'",  Becker*^,  Corum^,  Harris12®,  Klrstein1®®,  Singer2®®,  and 
q?8 

Weak-7  .  The  David  Taylor  Model  Basin  performed  several  tests  on 
controlled  spherical,  spec imens  1  ®® >  222 »  2®®  as  well  as 
cylinders12^2®®.  other  shapes  such  as  spheres1 ’ } 
ell  Ipso  Ids ^  9  >  33 ,  toroids1®®,  and  cones®0®-®0®  by  Holownia1**2 , 
Singer^1,  Williams-'1'7  and  Weller-7  have  been  studied  and  test 
data  ha-  been  compiled.  The  influence  of  stiffeners  has  also 
been  studied  experimentally.  The  influence  of  the  shape  and 
attachment  was  considered  by  Baruch®®,  Dow®®,  Card^2,^3,  Corum®®, 
Medgley20**  ,  S i nger‘‘®°  ’  2  ^®  and  Weller®20-®2®. 

Because  the  factors  that  influence  buckling  of  shells  have 
not  been  completely  understood,  experiments  have  been  performed 
and  data  collected  that  presents  conflicting  results.  The  col¬ 
lection  of  data  for  the  ratio  of  buckling  pressure  to 
classical  (P/P  )  for  the  clamped  spherical  cap1®1  having  various 
rise  heights,  \.  Figure  19  is  an  example  of  the  scatter  that  can 
be  anticipated.  The  size  of  the  specimen  is  usually  restricted 
by  the  load-carrying  capacity  (more  importantly,  the  stiffness) 
of  the  supporting  test  fixtures.  To  reduce  this  influence,  non- 
metallLc  test  specimens  have  been  used.  PVC,  acrylics,  rubber, 
mylar,  and  other  plastic  materials  have  been  candidates.  For 
thin  metal  shell  structures,  the  electrof orrning  process  has  also 
been  used76 , B58-P60. 
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The  unfortunate  thing  is  that  the  subscale  tests  are  also 
subject  to  the  same  physical  laws  as  their  prototypes.  This 
means  that  the  same  seven  control  parameters  suggested  by 
Thompson  apply  to  the  experimental  specimens  as  well.  Further, 
the  scaled  tests  may  amplify  one  imperfection  parameter  that  is 
not  actually  contributing  in  the  full-scale  specimens.  Since 
shape  imperfections  are  generally  used  as  one  of  the  controlling 
parameters,  special  care  should  be  taken  to  insure  that  only 
shape  imperfection  is  being  introduced.  The  control  of  these 
parametric  deviations  becomes  more  critical.  The  main  unknown 
factor  in  specimen  evaluation  is  the  determination  of  the  true 
state  of  stress.  Small  deviations  in  material  thickness  or 
residual  stresses  may  mask  the  desired  results.  Further,  when 
plastic  materials  are  employed,  the  uniformity  of  material 
"elastic"  coefficients  may  be  difficult  to  achieve  or  accurately 
assess . 

One  of  the  major  problems  facing  the  experimentalist  inves¬ 
tigating  the  buckling  of  shells  is  that  of  fabricating  the  "per¬ 
fect"  test  specimen.  The  fact  that  there  are  some  shell/load 
combi  nations  (especially  for  large  R/t  ratios)  that  are  very 
Imperfection  sensitive  is  borne  out  by  the  studies  of  Carlson, 
Sendelbeck  and  Hoff^  from  which  Figure  20  is  extracted.  One  can 
observe  the  significant  scatter  obtained  for  different  specimens 
and  the  authors'  achievement  of  better  correlation  between  exper¬ 
imental  and  theoretical  values  of  shell  buckling  pressure  as 
their  fabrication  processes  improved  for  the  later  specimens. 

They  used  an  electroforming  process  to  manufacture  their  test 
specimens,  which  were  essentially  complete  spherical  shells. 

For  large  radius  to  thickness  ratios  (R/t  >  1000)  and  near 
perfect  complete  spherical  shells,  Carlson,  Sendelbeck  and  Hoff 
observed  that  most  of  the  shell  surface  becomes  unstable  at  the 
same  pressure;  and  dimples,  covering  much  of  the  surface,  develop 
simultaneously.  However,  for  more  severe  flaws,  single  dimples 
form  and  appear  at  the  individual  flaw  locations.  It  has  been 
only  in  the  last  few  years  that  investigators  have  tried  to 
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ascertain  the  nature  of  imperfections  in  fabricated  shell  struc¬ 
tures,  either  full-scale  or  model-sized  (Arbocz^-^,  Horton^,,_ 

ancj  others » ^8  ) .  In  1969,  Arbocz  and  Babcock‘S  reported 
on  the  results  of  buckling  experiments  in  which,  for  the  first 
time,  the  actual  initial  imperfections  and  pre-buckling  growth  of 
electroplated,  isotropic  shells  were  measured  and  recorded  by 
means  of  an  automated  scanning  system.  In  1971,  Singer,  et  al., 
published  the  results  of  similar  imperfection  surveys  on  ring  and 
stringer-stiffened  shells.  Parmerter^®  in  _972  published 
results  on  an  interesting  holographic  method  of  recording  shell 
deformations  for  shallow  sperlcal  caps.  The  work  of  Krenzke  and 
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Klernan1  published  in  1965  addresses  the  problem  of  imperfec¬ 
tions  in  partial  nonshallow  spherical  shells  with  Imperfections 
machined  into  the  specimen. 

With  this  experience,  Horton^1*  ’  examined  the  collapse  of 
an  axially  loaded  cylinder.  It  was  suggested  that  as  the  load 
increased,  the  number  of  buckles  increased  with  the  change 
following  a  normal  Gaussian  distribution.  Thus,  Horton  concluded 
that  a  random  variation  In  the  imperfections  had  to  have  existed, 
[f  Thompson's  supposition  were  considered,  it  would  suggest  that 
the  number  of  equilibrium  paths  are  so  close  that  the  randomness 
would  be  attributed  to  the  number  of  possible  equilibrium  paths. 
Singer  has  shown  that  random  type  imperfections  can  be  correlated 
to  linear  buckling  theory.  Arbocz  has  demonstrated  that  with 
care  in  the  description  of  the  imperfection,  nonlinear  theory  can 
be  applied  to  give  very  accurate  results.  Therefore,  we  have 
contradicting  statements  for  the  same  phenomena. 

Babcock^  gives  a  rather  comprehensive  assessment  of  the 
experimental  methods  employed  from  1950  to  1972.  He  makes  a 
clarion  call  for  experimenters  to  carefully  design  the  test  and 
to  record  all  detailed  information  on  the  shell  specimen  and  its 
loading  environment.  The  call  appears  to  have  gone  unheeded. 
Because  the  archive  publication  requirements  have  restricted  the 
length  of  many  articles,  the  necessary  detailed  information  on 
specimen  configuration,  etc.,  must  usually  be  omitted  from  the 
article.  By  the  time  an  inquiring  analyst  pursues  the  detailed 


material,  the  data  has  been  mislaid,  or  the  data  is  part  of  a 
larger  set  of  data  which  is  under  a  corporate  proprietary  classi¬ 
fication  and  cannot  be  distributed,  or  the  researchers  have  left 
the  institution. 

Fortunately,  there  are  exceptions.  Arbocz  has  continued  his 
studies  of  Imperfection  in  cylindrical  shells  and  has  presented 
data  in  a  variety  of  publications  while  pursuing  imperfection 
correlation  studies.  Singer  and  Tennyson  and  their  associates 
also  have  contributed  to  the  traceable  effects  of  experimentally 
determined  imperfections  and  boundary  conditions. 

With  more  control  on  the  specimen  preparation,  it  is  possi¬ 
ble  to  examine  more  closely  the  effects  of  imperfection.  To  this 
end,  one  method  has  been  developed  to  prepare  specimens  with  a 
known  imperfection.  Boros^  examined  the  deteriorating  effect  of 
prescribed  axisymmetric  imperfection  shapes  on  the  buckling  load 
of  ring  and  stringer  stiffened  cylindrical  shells  under  a  variety 
of  loading  conditions.  He  developed  a  process  that  could  intro¬ 
duce  prescribed  axisymmetric  imperfections  into  stringer  stif¬ 
fened  cylindrical  shell  specimens  and  maintain  responses  in  the 
elastic  range.  Thus,  repeatability  of  results  could  be  demon¬ 
strated.  The  results  tend  to  show  that  Hutchinson's  nonlinear 
theory  is  quite  applicable  for  a  variety  of  shell  geometries  and 
stiffness  ranges  of  external  or  internal  stiffeners. 

In  a  recent  survey  on  buckling  experiments  of  shells,  essen¬ 
tially  conducted  after  Babcock's  review.  Singer1 2 3"^  reviewed  the 
relationship  between  the  need  for  experiment  and  the  use  of  the 
modern  digital  computer.  He  suggested  eight  reasons  why  experi¬ 
ments  of  the  buckling  shells  should  be  conducted: 

(1)  Better  understanding  of  buckling  behavior  and  the 
primary  factors  affecting  it. 

(2)  Finding  new  phenomena. 

(3)  Obtaining  better  inputs  for  computations. 

(A)  Obtaining  correlation  factors  between  analysis  and 
test,  and  for  material  effects. 
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(5)  Building  confidence  in  multi-purpose  computer 
programs. 

(6)  Testing  novel  ideas  of  shell  construction  or  very 
complicated  shell  elements. 

(7)  Understanding  buckling  under  dynamic  loading  and  in 
fluid/structure  interaction  problems. 

(8)  Certifying  tests  of  full  scale  shell  structures. 

With  this  survey  of  data  developed  in  the  70' s.  Singer  outlined 
eleven  areas  for  further  research,  including  many  of  the  areas 
covered  in  the  present  survey:  composites,  load  interaction, 
imperfection  data  banks,  mode  interaction,  boundary  conditions, 
and  testing  techniques.  It  is  recognized  that  although  adequate 
computational  methods  are  available  for  the  theory  upon  which 
they  are  based,  a  need  for  the  verification  of  the  theory  itself 
remains  a  requirement.  This  can  only  be  accomplished  through 
verification  by  experiment. 
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6.0  CONCLUSIONS  AND  SUMMARY 

It  appears  that  the  state-of-the-art  of  the  analysis  of 
imperfection  sensitivity  of  shells  is  not  at  a  sufficient  level 
to  provide  an  all-encompassing  theory  that  engineers  can  use  with 
absolute  confidence.  The  theory  of  Koiter  has  been  generally 
accepted  as  one  measure  in  determining  the  stability  character¬ 
istics  of  shells  but  this  theory  has  proven  to  be  limited. 
Extensions  of  this  theory  by  Thompson,  et  al. ,  have  shown  some 
promise  to  extend  the  range  of  shell  properties  over  which  their 
method  could  be  applied,  but  more  importantly  they  make  it  pos¬ 
sible  to  identify  buckling  phenomena  that  could  be  overlooked  by 
others. 

As  pointed  out  in  this  survey,  the  theorectical  development 
on  the  understanding  of  imperfection  sensitivity  of  shells  can 
become  quite  cumbersome  even  when  examining  the  behavior  of 
simple  shapes.  The  developments  of  the  Koiter  postbuckling  pro¬ 
cedure  falls  into  this  category.  The  conditions  set  forth  by 
Thompson  is  almost  of  the  same  complexity.  However,  there  seems 
to  be  some  benefit  in  developing  an  understanding  for  both 
methods.  If  semi-empirical  methods  are  to  be  used,  there  is 
always  a  question  as  to  the  condition  of  the  shell  under  investi¬ 
gation  as  compared  to  the  accepted  test  data.  Further,  one  must 
accept  the  lower  bound  philosophy  inherent  with  semi-empirical 
methods . 

As  an  aid  to  the  reader,  a  partial  categorized  cross- 
reference  is  provided  in  Table  1  for  the  various  articles 
highlighted  in  this  survey.  Their  order  of  appearance  by  no 
means  implies  a  rating  of  the  material.  The  categories  are 
subdivided  as  general  theory  of  shells,  state  of  art  review, 
stability  prediction  methods,  numerical  methods  used  in  shell 
stability  analysis  and  experimental  results  of  different  shell 
geometric. 

At  present,  only  one  computerized  procedure  (written  by 
Cohen)  permits  the  study  of  imperfection  sensitivity  of  an  arbi¬ 
trary  shell  of  revolution  under  specific  loads  using  Koiter* s 
theory.  In  order  to  perform  a  simple  investigation,  Cohen's 


TABLE  1.  CROSS-REFERENCES  BY  CATEGORY 


GENERAL  THEORY 


TABLE  1.  (continued) 


NUMERICAL  METHODS 


Finite  Element 

Finite  Difference 

Numerical  Integration 

122-123,  199 

7,  43,  66,  80 

165 

83-86,  87 

Experimental 


Cylinder 

2,  4,  12-18,  23 

30, 

36,  46,  53 

54, 

69,  72-75 

88, 

90,  100,  128 

129 

132,  139 

140 

162',  196 

204, 

213-215 

231- 

-236,  238 

244, 

249-253 

259, 

282,  290 

305, 

320-327 

328 

130,  142 
145,  260 
263,  267 
329 


Spherical 


76,  159,  167 
183,  184,  186 
193,  202,  222 
256 


19,  33,  77 
111,  141,  15 
186,  296 


procedure  requires  the  use  of  three  separate  computer  programs. 
This  approach  has  not  been  widely  used  because  of  its  complexity 
and  availability.  Other  computer  programs  that  can  compute  the 
buckling  load  of  imperfect  shells  are  available,  but  the  cost  of 
the  analysis  may  become  prohibitive. 

It  is  comforting  to  note  that  when  these  numerical 
procedures  are  applied  to  specific  controlled  problems,  limited 
correlation  between  the  previously  mentioned  theories  and  the 
more  general  computerized  procedures  have  been  made.  If  the 
analyst  is  willing  to  accept  the  results  from  these  procedures, 
he  must  also  be  aware  of  their  shortcomings.  The  stability  char¬ 
acteristics  of  the  structure  are  being  evaluated,  not  the  numer¬ 
ical  procedure.  Further,  since  most  computer  procedures  have 
certain  computational  capacity  limitations,  the  principal  compo¬ 
nents  that  will  develop  critical  participating  responses  must  be 
approximately  known  in  advance.  Thus,  non-critical  components 
need  not  be  considered  in  the  formulation  of  the  problem  but 
there  is  always  a  chance  of  not  including  a  principal  component. 

The  need  for  a  simple  computer  program  for  the  simple  evalu¬ 
ation  of  imperfection  sensitivity  is  required.  This  can  be  done 
for  uniform  shells,  such  as  cylinders,  but  not  for  general  shells 
or  loads.  Further,  the  study  of  stiffened  shells  with  imperfec¬ 
tions,  which  can  be  performed  on  a  simplified  basis,  needs  to  be 
made  available  to  the  engineer.  With  these  tools,  the  engineer 
can  at  least  make  some  preliminary  evaluations  without  the  need 
of  complex  relationships  or  large  computer  budgets  being 
considered. 

More  general  or  unified  theories  are  required  to  handle  the 
specific  case  studies  required  in  the  analysis  of  near  optimum 
shell  structures.  This  can  only  be  accomplished  through  the 
extension  of  the  accepted  theories  and  the  accompanying  numerical 
methods.  Again,  as  pointed  out  by  Singer,  improvement  in  the 
accuracy  of  these  prediction  techniques  is  required,  and  they 
must  be  supplemented  and  verified  through  adequate  experiments. 
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DISCUSSION 

■I.  Singer,  Techlon  Israel  Institute  of  Technology 

The  review  is  a  thorough  and  comprehensive  survey  of  the 
field,  attempting  to  bridge  the  gap  between  basic  shell  theories 
and  practical  elastic  stability  analysis  for  stiffened  shells. 
This  difficult  task  is  accomplished  successfully  in  most  of  the 
review,  but  the  following  comments  are  warranted: 

1.  The  possible  practical  applications  of  catastrophe 
theory  are  not  obvious  from  the  brief  discussion  presented,  and 
the  implications  of  the  two  theorems  require  further  clarifica¬ 
tion. 

2.  The  effect  of  boundary  conditions  on  unstiffened  cylin¬ 
drical  shells  is  briefly  mentioned  in  the  review.  For  stiffened 
shells,  however,  the  influence  of  the  boundary  conditions  is  of 
prime  importance  and  may  sometimes  even  overshadow  the  influence 
of  geometrical  imperfections  (see,  for  example.  Refs.  273,276,  or 
Refs.  A1-A3*),  and  hence  justify  some  elaboration.  The  influence 
of  load  eccentricity  may  also  be  of  prime  importance  (see  Refs. 

51  or  324). 

3.  The  predominant  effect  of  the  boundary  conditions  on  the 
buckling  of  stiffened  shells,  and  their  similar  influence  on  the 
vibrations  of  these  shells,  motivated  extensive  correlation 
studios  which  yielded  a  nondestructive  tool  for  definition  of  the 
boundary  conditions  —  the  vibration  correlation  technique.  This 
technique,  which  consists  essentially  of  an  experimental 
ieterminat ion  of  the  lower  natural  frequencies  of  a  loaded  shell 
and  evaluation  of  equivalent  elastic  restraints  that  represent 
the  actual  boundary  conditions,  has  been  developed  in  recent 
years  and  applied  to  different  shells  and  loading  conditions  (see 
Refs.  A4-A8 ) . 


#The  prefix  "A”  refers  to  Additional  References  following 
Author's  Closure. 
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The  imperfection-sensitivity  predictions  shown  in  Pigs. 
16-18  of  the  review  are  found  in  later  studies  (see  Ref.  153)  to 
depend  strongly  on  the  prebuckling  deformations,  which  can 
appreciably  change  the  imperfection-sensitivity  and  the  7. 

(Batdorf  shell  geometry  parameter)  at  which  it  occurs  (see 
Fig.  11  in  Ref.  268).  Furthermore,  tests  do  not  verify  the 
predictions  of  Figs.  16-18  (see  Fig.  14  of  Ref.  A3). 

5.  Extensive  imperfection  measurements  have  recently  been 
carried  out  on  stiffened  shells  and  an  international  data  bank  of 
imperfection  measurements  has  been  established  (see  Refs.  26,28 
and  A9 )  which  will  eventually  permit  meaningful  correlation 
studies  between  manufacturing  techniques  and  significant 
Imperfection  shapes. 


A.  van  der  Neut,  Delft  University  of  Technology 

The  survey  deals  with  an  enormous  amount  of  literature  and 
the  author  has  done  a  very  good  job.  My  response  on  the  parts 
with  which  the  writer  is  familiar  are  noted. 

The  usual  approach  of  establishing  the  minimal  equilibrium 
load  at  large  deformations  was  known  before  the  war,  due  to  von 
Karman  and  Tsien's  paper  of  1939  (see  Reference  314).  Therefore, 
Hotter  deliberately  rejected  this  approach  when  he  chose  to  focus 
on  the  circumstances  which  trigger  collapse  and  on  the  load  at 
which  collapse  occurs.  The  difference  between  the  two  philoso¬ 
phies  can  be  illustrated  by  the  inves tigation  of  a  car  accident 
where  the  car  crashed  into  a  tree.  The  parallel  to  establishing 
the  far  post-buckling  equilibrium  is  then  the  assessment  of  the 
amount  of  damage  done  by  the  tree,  whereas  the  real  problem  is  to 
know  what  caused  the  car  to  deviate  from  its  straight  course. 
During  two  decades  and  even  during  the  sixties,  when  Koiter's 
theory  had  received  publicity,  a  vast  amount  of  energy  was  wasted 
In  studying  the  final  "car  damage"  before  it  was  generally  recog¬ 
nized  that  the  real  problem  is  the  initial  post-buckling 
behaviour. 
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The  two-f lange-model  the  writer  used  in  Ref.  310  was 
part  1  cul  nrly  apt  to  call  attention  to  the  unfa  vourab  Le  effect  of 
interaction  between  overall  and  local  buckling,  because  it 
exaggerates  this  effect  due  to  the  presence  of  two  independent 
local  buckling  modes  at  the  same  stress.  A  Later  report  (see 
Reference  A10)  suffers  from  the  same  defect. 

An  alternative  approach  to  the  interaction  between  local  and 
overall  buckling  in  stiffened  panels,  given  by  W.T.  Koiter  and  M. 
Tignataro,  avoids  the  coincidence  of  two  local  modes  and  yields  a 
milder  effect  of  interaction.  In  later  work  the  writer  removed 
the  defect  of  coinciding  local  modes  by  introducing  elastic 
coupling  between  the  two  flanges.  It  lessened  the  effect  of 
interaction  and  confirmed  Koiter* s  "alternative  approach"  for  the 
case  where  the  stiffener  walls  in  local  buckling  have  much 
smaller  deflections  than  the  skin  plate.  It  appeared  however 
when  applying  the  "alternative  approach"  and  the  revised  model  to 
a  panel  where  the  deflections  of  the  stiffener  top  (e.g., 

Z-st  if  Zoned  panel)  are  of  the  same  order  of  magnitude  as  those  of 
the  skin ,  the  Koiter  method  underestimates  the  effect  of  inter¬ 
action.  Thereupon,  Koiter  has  Improved  his  analysis  by  taking 
into  account  the  first  two  local  modes  instead  of  only  the  first 
one.  The  revised  model  parameters  were  adjusted  to  the  results 
obtained  by  Koiter  for  the  perfect  structure.  This  model  was 
then  used  to  establish  the  effect  of  local  imperfections.  Next, 
for  the  square  tube  the  stability  at  the  bifurcation  point 
against  small  finite  deflections  of  the  column  axis  and  the 
effect  of  columns  axis  imperfections  were  established.  These 
most  recent  developments  have  been  the  subject  of  a  joint  paper 
by  Koiter  and  the  writer  and  were  presented  at  the  "Int.  Conf.  on 
Thin-Wal  led  Structures"*'®'*'. 


AUTHOR'S  CLOSURE 

I  wish  to  thank  Profs.  Singer  and  van  der  Neut  for  their 
added  discussions  on  the  survey.  Since  a  great  deal  of  informa¬ 
tion  had  to  be  omitted  from  this  survey,  some  of  the  key  points 
that  describe  the  basic  assumption  behind  the  theories  may  not  be 
sufficiently  highlighted.  To  the  reader,  I  apologize. 

Singer's  first  remark  with  respect  to  the  use  of  catastrophe 
theory  should  be  addressed.  This  very  point  has  raised  a  few 
negative  remarks  by  other  researchers. 

Thompson  and  Hunt  in  their  book  present  a  rather  detailed 
discussion  of  how  energy  methods  (not  catastrophe  theory)  can  be 
employed  in  the  evaluation  of  elastic  stability.  The  two  stated 
theorems  permit  elastic  stability  prediction  to  be  made  solely 
from  equilibrium  studies.  Further  correlation  between  Koiter's 
theory  (using  energy  theories)  and  catastrophe  theory  (using 
topological  similarities)  has  been  suggested  by  Thompson,  (Ref. 
All).  Hanson,  (Ref.  A12),  has  demonstrated  that  through  an 
extension  of  Koiter's  method,  a  two-mode  buckling  problem  is 
shown  to  be  represented  by  a  three-dimensional  surface.  He 
classifies  this  surface  as  to  bifurcation  and  correlates  these 
with  catastrophe  theory.  This  is  again  illustrated  by  Hui  and 
Hansen,  (Ref.  A13),  in  a  recent  publication.  In  a  separate 
publication  Hansen  and  Hui  (Ref.  A14)  show  that  the  potential 
function  for  a  two  mode  buckling  of  a  spherical  shell  can  be 
described  as 

p  p  ii  p  p 

R  =  x  y  +  y  +  tx^  +  wyc  -  ux  -  vy 

where  x  and  y  are  amplitudes  of  certain  deflection  modes  and  u, 
v,  w  and  t  are  either  applied  loads  or  initial  imperfections. 

This  is  called  the  parabolic  umbilic  catastrophe.  By  applying 
the  theory  of  Thom  (Ref.  A15),  the  question  of  stability  can  be 
resolved.  Hansen  and  Hui  show  that  by  examining  the  stable  and 
instable  forms  of  the  above  function,  the  complete  behavior  of  a 
sphere  subjected  to  uniform  and  varying  pressures  can  be  mapped. 
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Alesso  [Al6]  has  Further  demonstrated  how  the  mapping  of  a 
^  |) 

control.  parameter  ("gy)  onto  a  three-dime;  tonal  surface  relates 
to  the  possible  deformations  that  can  occur  in  a  shallow 
spherical  cap.  Figure  A1  illustrates  the  interrelationship 
between  elementary  catastrophe  theory  and  that  suggested  by  the 
more  accepted  energy  methods.  Clearly,  the  two-dimensional  plane 
suggested  by  energy  methods  contains  less  information  than  that 
developed  from  catastrophe  theory. 

If  one  were  to  fully  understand  the  implications  of  the 
catastrophe  theory,  the  mathematical  description  of  the  potential 
function  developed  for  structural  stability  may  be  more  readily 
understood,  particularly  for  coincident  bifurcation  loads.  For 
single  buckling  paths,  catastrophe  theory  does  not  appear  to  have 
an  advantage  over  energy  methods.  The  major  advantage  that  can 
be  utilized  by  employing  catastrophe  theory  is  that  a  qualitative 
measure  of  the  effects  of  imperfections  can  be  obtained.  Which 
imperfection  shape  and  the  number  of  terms  required  for  proper 
convergence  can  be  identified  using  this  theory.  The  method  is 
not  quantitative.  The  key  here  is  understanding  some  new 
mathematical  development.  At  present,  the  newly  trained  design 
engineer  would  probably  not  have  the  discipline  nor  interest. 
However,  It  may  be  time  for  the  engineer  to  open  his  eyes  and 
possibly  develop  the  "by  hook  or  by  crook"  method  Koiter  was 
alluding  to. 
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